arXiv:hep-ph/0112035vl 3 Dec 2001 


The Top Quark at Linear Colliders: 
Quantum Chromodynamics Corrections 


by 


Cosmin Macesanu 


Submitted in Partial Fulfillment 
of the 

Requirements for the Degree 


Doctor of Philosophy 

Supervised by 
Professor Lynne H. Orr 

Department of Physics and Astronomy 
The College 
Arts and Sciences 


University of Rochester 


Rochester, New York 


2001 


11 


to my family 



Ill 


Curriculum Vitae 


The author was born in Bucharest, Romania, on May 13 1973. He attended 
the Faculty of Physics at the University of Bucharest from 1991 to 1996, and 
graduated with a Bachelor of Sciences degree in 1996. He came to the University 
of Rochester in the fall of 1996 and began graduate studies in the Department of 
Physics and Astronomy. He received a Robert Marshak Fellowship in 1996 and 
1997. He pursued his research in the field of elementary particle physics under 
the direction of Professor Lynne H. Orr and received the Master of Science degree 
in 1997. 



Acknowledgments 


First, I want to thank my adviser, Professor Lynne H. Orr, for her support 
of my work on this thesis and for suggesting the topic in the first place. 1 would 
also like to express my gratitude for the cooperation and help of my colleagues 
and of many postdoctoral fellows who performed research at the University of 
Rochester these past five years. Among these, I would like to give special thanks 
to Dr. Doreen Wackeroth, who spend many hours discussing with me the details 
of the evaluation of the electroweak corrections to the W pair production and 
decay process, and to Dr. Rob Szalapsky, who helped me with the computation 
of Passarino-Veltman functions and with the computational details of my work 
in general. Moreover, 1 would like to thank Dr. Oleg Yakovlev for a discussion 
concerning the diagrams contributing to the top production and decay process, 
and Dr. Alexander Chapovsky for discussions on the DPA approximation and for 
providing me with routines for the evaluation of nonfactorizable corrections in the 
soft gluon approximation. 



V 


The Top Quark at Linear Colliders: 
Quantum Chromodynamics Corrections 

by 

Cosmin Macesanu 


Abstract 

We present a computation of QCD next-to-leading order corrections to the top 
production and decay process at linear colliders. The top quarks are allowed to 
be off-shell and the production and decay subprocesses are treated together, thus 
allowing for interference effects. We consider the case of real gluon radiation, as 
well as virtual corrections to the tree level amplitude. The framework employed 
for our computation is the double pole approximation (DPA). We describe the 
implementation of this approximation for the top production and decay process 
and compare it with the implementation of DPA for the evaluation of QED cor¬ 
rections to the W pair production at LEP II. Similarities and differences between 
the two cases are pointed out. 

The theoretical approach we present in this thesis is implemented in a Monte 
Carlo generator. The total amplitude is separated into gauge invariant parts which 
can be associated with radiative corrections to production and decay subprocesses, 
and interference between these. The contributing amplitudes are computed using 



spinor techniques, and include all top width effects, spin correlations and b quark 
mass effects. The results discussed for the real gluon radiation case include studies 
of the gluon radiation properties and the effects of this radiation on top mass 
reconstruction. We also examine the effects of interference between production- 
and decay- stage radiation, whose magnitude is sensitive to the value of the top 
quark width. After the computation of virtual corrections, we present results 
for the total top production cross sections, and we analyze the magnitude of 
nonfactorizable (interference) corrections. We study the impact these corrections 
can have on the top invariant mass distributions. 
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Chapter 1 
Introduction 


Since ancient times, people have been looking for explanations for the natural 
phenomena surrounding us. With the birth of the modern scientific method in 
the 15th - 16th century, great steps forward have been made in our understanding 
of the natural laws. The body of accumulated knowledge has naturally coalesced 
into theories which deal with different aspects of reality. Examples of such theories 
are Newtonian mechanics, which deals with the interaction of normal bodies, 
optics, which deals with properties of light, and electricity and magnetism, which 
deals with the electric and magnetic properties of materials. As time has passed, 
there has been a tendency to look for ways to combine these disparate theories 
into more fundamental ones, which encompass and explain all the phenomena 
previously dealt with separately (Maxwell’s combination of electricity, magnetism 
and optics into electromagnetism is such an example). In the past century, these 
searches have coalesced into a search for a unified theory which underlies all the 
physical reality around us. 

The result is the Standard Model of particle physics. This theory deals with 
the fundamental constituents of matter and their interactions (except gravity). 


1 
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According to our present knowledge, all ordinary matter is built from elementary 
particles: six leptons and six quarks (spin 1/2 particles, or fermions), which are 
arranged into three families: 


Leptons : 


Quarks : 



(i.i) 


The interactions of these particles are mediated by the gauge bosons: 


Gauge bosons : 7 , Z 0 , W ± , g 


( 1 . 2 ) 


which are spin 1 particles. The mathematical framework which describes these 
interactions is gauge field theory; according to this theory, each generator of the 
gauge symmetry group of the Lagrangian corresponds to one gauge force carrier. 
Thus, the photon ( 7 ), Z 0 and are carriers of the electroweak gauge force 
(symmetry group U( 1) x SU(2)), which is felt by all particles, while the gluons 
(g) are the carriers of the strong gauge force (symmetry group SU(3)) which 
mediates the interactions of the quarks among themselves solely. 

So far, the successes of the Standard Model are impressive; some of its predic¬ 
tions have been tested with a remarkable degree of accuracy, and, at this time, 
not a single piece of experimental evidence contradicts it convincingly. Still, there 
are pieces missing; for example, the exact mechanism of electroweak symmetry 
breaking (which is responsible for giving masses to the particles in Eqs. |1 . 1| , 0 is 
still unknown. It is surmised that the breaking of electroweak symmetry is driven 
by yet undiscovered bosons: either the Higgs particle(s), in the most promising 
models, or maybe top quark condensates (in technicolor models [jlj). Supersym¬ 
metric (SUSY) models 0, which solve some theoretical problems in the Standard 
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Model, also predict a whole slew of new particles (supersymmetric partners). In 
general, there is a great interest today in what lies beyond the Standard Model. 
(Since it does not include gravity, we know that it cannot be the final theory). It 
is expected that, in coming years, experiments at higher energies will help answer 
these questions. 

The study of the top quark might shed light on the answer to at least some 
of these questions. The characteristics of the top quark make it one of the most 
interesting elementary particles discovered so far. Its mass is quite large, about 
175 GeV (we use natural units in which h = c = 1); correspondingly, its Yukawa 
coupling (the coupling to the Higgs boson) is of order unity; and this might be an 
indication that the top plays a special role in the clectroweak symmetry breaking 
process. Moreover, the top Standard Model width is about 1.5 GeV, thus being 
much larger than A qqd ~ 200 MeV. This means that the top quark decays before 
having time to hadronize ||, therefore providing us with a unique opportunity 
to study the interactions of a bare quark. All these properties, together with the 
fact that, since the energies involved in the top production and decay processes 
are large, we can use perturbative QCD for reliable theoretical predictions, insure 
that the study of the top quark will be one of the main goals of particle physics 
for the next decade. 

The top quark was discovered in 1995 at the Fermilab Tevatron Collider [Q. 
Only a couple hundreds events have been identified so far. Due to limited statis¬ 
tics, the only precise information available on the top quark so far is its mass Q; 
the latest analysis || gives the value 174.3 ±5.1 GeV. The Tevatron Run II, start¬ 
ing this summer, is expected to provide us with a sample of top events about an 
order of magnitude larger than the one available so far; this will result in more 
1 There are also results for other top quark parameters, like the production cross section and 
its couplings, but the values obtained have large statistical uncertainities. 
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precise determination of the top mass, width and couplings. The Large Hadron 
Collider, once it starts operating in 2006, will be a top factory, producing more 
than 8 million top-antitop pairs anually 0. However, the analysis of the data 
coming from hadron colliders is complicated by uncertainties in the initial state, 
large QCD backgrounds, etc. Thus, even with the large top sample provided by 
the LHC, the uncertainity in the top mass measurement, for example, will be of 
order 1 to 2 GeV (due mostly to systematic uncertainties 0). 

An e + e _ linear collider with center of mass energy greater than the tt produc¬ 
tion threshold (350 GeV) would be an ideal machine for the precision study of 
the top quark (among other things). At such a collider, the electron and positron 
in the initial state annihilate and create a top-antitop pair, which in turn each 
decay into a Wb pair. So far, this machine is in the design stage, with a projected 
date of completion not before the year 2010. At this time, there are two main 
designs; the European TESLA and the American-Japanese NLC-JLC (for details 
about the parameters of these machines, see for example 0). Both machines are 
linear accelerators, about 30 Km in length, operating at energies up to 1 TeV 
(800 GeV for TESLA) with luminosities of order 10 34 cm _1 s _1 (which means 100 
inverse femtobarns per year). 

One of the main advantages of an e + e~ collider versus a hadronic one is that 
the energy of the initial state is well determined. Thus, it is possible to perform 
threshold studies for the production of a tt pair. The shape of the production cross 
section near the threshold is quite sensitive to a number of top quark parameters. 
In this energy range, current analysis indicates that it is possible to measure the 
top mass with an accuracy of about 40 MeV using only 10 fb~ l 0. With a larger 
data sample, measurements of the top decay width, Yukawa coupling and strong 
coupling constant at the several percent level can also be achieved. Going to 
higher energies, we can study the V-A structure of the top quark couplings to the 
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gauge bosons ( 7 , Z and W ) The information about couplings can be obtained 

by using spin correlations: the top quarks are produced in certain spin states, as 
dictated by the top - 7 , Z couplings. Since the top decays before hadronization, 
the spin states of the top directly influence the angular distributions of its decay 
products. Simulations show that by analyzing kinematic variables of final state 
particles we can measure top couplings at the several percent level 0 . 

Of course, in order to obtain information about fundamental parameters like 
the top quark mass and couplings, experimental data is only a part of the equa¬ 
tion. The other part is a good theoretical understanding of the physical processes 
which are studied. For the threshold region, comprehensive theoretical studies 
(NNLO computation with resummation of large logarithms, careful treatment of 
the renormalon ambiguity) have already been performed 0 . Above threshold, the 
theory lags behind. At 500 GeV center of mass energy, with a 500 pb~ l integrated 
luminosity, we will get about 3 x 10 5 top quark pairs created; as a consequence, 
the experimental accuracy is better than one percent. Ideally, we would like to 
have a similar (or better) precision in the theoretical predictions. 

As for most physical processes, precise theoretical predictions for the top pro¬ 
duction and decay require the computation of next to leading order corrections 
(either virtual or real). These can be split into two categories: electroweak cor¬ 
rections, due to the radiation of a photon or weak interaction gauge boson (Z or 
IT), and QCD corrections, due to virtual or real gluon radiation. Since the top is 
a quark, the QCD corrections are the most important ones, and these corrections 
will be the subject of the present paper. 

First, we start by reviewing previous work done in the area. In the first ap¬ 
proximation, the top production and decay processes can be considered separately. 
The QCD radiative corrections to individual subprocesses have been comprehen- 
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sively studied. For the top production subprocess: 


e e 


t t C 9 ) 


there are computations for virtual and soft (low energy) gluon radiation m , as 
well as for hard (higher energy) gluon radiation ( 0], 0, B 0 are just some 
examples). Similarly, the top decay subprocess: 

t^bW(g) 

has been computed taking into consideration virtual and hard gluons together 

(0,0, H). 


Using these results, one can try to approximate the top production and de¬ 
cay process by assuming that the intermediate tops are on the mass-shell (narrow 
width approximation) and treating the subprocesses separately Jl7| . This assump¬ 
tion is usually reasonable; the result for the total cross section is valid up to terms 
of order T t /m t ~ 1%. However, finite top width effects, which can be thought of 
as interference between production and decay subprocesses, can be important in 
some differential cross sections. Therefore, for precision better than % level, it is 
necessary to treat the production and decay together, by allowing the top to go 
off-shell. Thus, in our computation we take into account the interference terms 
(also known as nonfactorizable corrections). 

We shall perform this computation and present the results at the parton level 
only. Since the final state of the top quark pair production process is bW + bW~ 
(at lowest order), the experimental signature is either 6 jets (if both W’s de¬ 
cay hadronically) or lepton(s) + jets + missing energy (if one or both W’s have 
semileptonic decays). We assume that the issues related to jet reconstruction and 
identification have been solved, and our final state contains two W bosons fj. two 
2 At the Monte Carlo level, we actually allow the on-shell W’s to decay, either semileptonically 
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b quarks and possibly a gluon. Even at this level, the complete computation of 
all the diagrams contributing to this final state (Born and next to leading order) 
is a very difficult task. Therefore, we shall employ the double pole approximation 
(DPA) which means taking into account only the diagrams which contain a top - 
antitop pair. 

It is worth mentioning a resemblance between the process of interest to us: 
e+e" —>tt(g)^b W + b W~ (g) 


and the W pair production and decay process (e + e~ —> W + W~ —► 4/) with 
QED corrections at LEP. The issues which arise in the two computations are 
similar, because in both cases we are dealing with the production and decay of 
heavy unstable particles. Our treatment is largely similar to the one used for 

m. B- But there are some differences, both in 


the electroweak process 


the implementation of the DPA approximation and in the number and type of 
terms which contribute to the final result (the latter being due to the fact that in 
our case the intermediate off-shell particles are fermions, and not bosons). These 
differences will be pointed out in the course of our discussion. 

The outline of the thesis is as follows. In Chapter 2 we lay out the general 
framework in which we perform our computation. This includes a description of 
the double pole approximation (DPA) method for the top production and decay 
case, with an overview of the Feynman diagrams contributing to this process. We 
also discuss the computation of differential cross sections, with some details on 
the treatment of the infrared singularities, and we review some salient features of 
the treatment of the widths of unstable particles. 

Chapter 3 is dedicated to an analysis of the top production and decay process 


or into a pair of massless quarks, but in the latter case we do not take into consideration QCD 
corrections to the W decay. 
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in conjunction with the radiation of a real (hard) gluon. We start with a de¬ 
scription of the computation of the amplitudes contributing to this process. Since 
we perform a full off-shell computation, a prescription for separating the total 
amplitude into parts which can be associated with radiation in the production or 
decay stage is also given. Finally, we study the properties of gluon radiation, and 
the impact of the gluon on top mass reconstruction. Special attention is paid to 
the issue of interference between production- and decay-stage radiation, and the 
effect which it can have on differential distributions. 

The computation of virtual corrections is the subject of Chapter 4. The eval¬ 
uation of NLO amplitudes in the double pole approximation is presented and 
discussed in some detail. The results for the amplitudes corresponding to interfer¬ 
ence diagrams are similar to results previously obtained for the W pair production 
process, while for the vertex corrections and fermion self-energy diagrams there 
are differences between the two cases. In order to facilitate comparisons with 
the on-shell approach, we also formulate our results in terms of correction to the 
production and decay subprocesses and interference contributions. The gauge in¬ 
variance of the total and partial amplitudes in DPA is manifest in this formulation. 
Results for this section include values for the total cross section and an analysis 
of the impact of interference effects on the invariant top mass distribution. Com¬ 
parisons between our results and results obtained in an alternative approach (in 
which the real gluon interference terms are computed analytically) are also made. 

We end in Chapter 5 with a summary and the conclusions. The Appendix 
contains some technical details related to our computation: we present a short 
description of the spinor techniques used to evaluate the amplitudes contributing 
to our process, expressions for these amplitudes, and details about the evaluation 
of virtual corrections with the help of Passarino-Veltman functions. 



Chapter 2 

Computational Approach 


Our aim in this thesis is the evaluation of next-to-leading order amplitudes and 
differential cross sections for the top production and decay process at linear col¬ 
liders. In this section, we start with a description of the general framework in 
which we perform our computation. The Feynman diagrams contributing to the 
process are introduced, and the double pole approximation is explained. Some 
general issues related to the evaluation of the cross section, the treatment of in¬ 
frared singularities and the treatment of the top width are examined. Details 
about the evaluation of amplitudes and a comprehensive analysis of the specific 
issues arising in this computation can be found in later sections. 

2.1 Amplitudes in DPA 

In top pair production at linear colliders, what is actually observed experimentally 
is the process 

e+e" -> bW + bW~ (2.1) 


9 



CHAPTER 2. COMPUTATIONAL APPROACH 


10 


There are many diagrams contributing to this process. At tree level, they can be 
split into 3 classes: diagrams which contain a top-antitop pair (Fig. |2-1|) , dia¬ 
grams which contain either one top quark or a top antiquark (Fig. \2-2\ and charge 
conjugates), and diagrams which do not contain any top (there are about 50 such 
diagrams). Although not an easy task, it is possible to perform the evaluation of 
all these amplitudes (by using one of the automated tree level amplitude compu¬ 
tation programs, like MADGRAPH |22|]). The computation of QCD corrections 
to all tree level diagrams increases the degree of complexity by quite a lot, and is 
probably not feasible yet. 



Figure 2-1: The top-antitop diagrams contributing to the process (|2 . 1|) . 


Fortunately, we can make use of the fact that, looking for top pair production, 
we are interested in a specific region of the final state phase space of ( |2.1|) . This 
region is defined by the requirement that the invariant mass of the W b pair is close 
to the top mass: pj , p f ~ rnf , where p t = Pw+ + Pb and pi = Pw~ + Pi D I n this 
region, the amplitudes corresponding to the top-antitop diagrams are enhanced 

1 pi are the four-momenta of the particles involved in the process; = (E-,p x ,p y ,p z ) and 
p 1 2 = p^Pfi = E 2 — p x — Py — pi is a quantity invariant under Lorenz transformations. We also 

denote the antiquarks by a bar over the corresponding quark symbol. 
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Figure 2-2: Single top diagrams contributing to the process fl2Tj) . 


by the two resonant propagators coming from the two intermediate top quarks: 

1 1 


M 


2 — 2 2 — 2 
PI ~ mi pf - mi 


( 2 . 2 ) 


For this reason, the first class of diagrams (Fig. E3) are called doubly resonant 


diagrams. Correspondingly, the second class (Fig. |2-2|) and the third class of 
diagrams, which contain only one top quark propagator or none, are called singly 
resonant and non-resonant diagrams respectively. 

Note that the contributions coming from singly resonant diagrams (whose 
corresponding amplitudes contain a single resonant propagator) are reduced by a 
factor T t /m t with respect to the doubly resonant contributions. Therefore, in the 
first approximation we can neglect the singly-resonant and non-resonant diagrams; 
we shall use the double pole approximation (DPA) which means keeping only the 
amplitudes which have a doubly resonant behavior when the top and the antitop 
go on-shell. 

At next to leading order, the diagrams contributing to the top production and 
decay process are presented in Figures |2-3| and [2-4| . In Figure [2-3| we present the 
doubly-resonant diagrams contributing to the process with a real gluon in the final 
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b 


b 




Figure 2-3: Feynman diagrams for gluon radiation in top production and decay. 

state: 

e+e^ -> bW + bW~g (2.3) 

The DPA amplitude for this process can be written as a sum of four terms : 

M ra = M t + Mi + M h + M- b (2.4) 

each term corresponding to one of the diagrams in Figure [2-3| (the subscript refers 
to the quark the gluon couples to). 

Some of Feynman diagrams for virtual corrections to the tree level process are 
presented in Figure [2A| . These diagrams can be divided into two classes : correc- 
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tions to particular subprocesses - the vertex and fermion self-energy diagrams in 
Fig. |2-4| a) and b) respectively - and interference type corrections (Fig. |2-4| c) 
and d)). Strictly speaking, the vertex and self-energy diagrams also contribute to 
interference between subprocesses; but, for the sake of brevity, we shall refer to 
the diagram in Fig. |2-4j a) as the production vertex correction diagram, and so 
on. Also, in the following, we will denote the tree level amplitude (Fig. |2-1|) by 
AA°, and the amplitude for the first order virtual corrections by A4 V9 : 

M vg = A4tt + AAtb + A4^ + AA.u + AA tb + A4 bb (2-5) 

ffere, the first three terms correspond to the three off-shell vertex corrections 
(which include in a suitable way the fermion self-energies, as described in section 
4.1.3), and the last three terms come from the interference diagrams. 

2.2 Cross sections and infrared singularities 

The partial amplitudes corresponding to the contributing Feynman diagrams are 
evaluated using spinor techniques (see Appendix A for details). The differential 
cross sections are obtained by taking the square of the total amplitude (sum of 
partial amplitudes) and summing over the polarizations (spins) of particles in the 
final state: 

da°’ V9 ’ r9 oc ^ |M°™| 2 (2.6) 

polarizations 

The proportionality constant contains phase space and kinematic factors and it is 
given in Appendix B. Our approach allows the computation of cross sections for 
particles of definite polarizations in the final state; but since these polarizations 
cannot be experimentally observed, all our results will be given with the spin sums 
in the final states performed. 
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Figure 2-4: Feynman diagrams for virtual gluon corrections to top production and 
decay. 

Our computations are relevant to two different experimental situations: the 
case when the final state contains only the b, b, W + and W~ particles, and the 
case when besides these particles there is also a gluon. In the latter case, the 
only amplitudes which contribute are the amplitudes in Figure |2-3|; the total cross 
section can be written as: 

a r9 cx [ \M r9 \ 2 dO 0+g (2.7) 

(where h2 0 is the phase space of the tree level process, and h2 0 +g is the phase space 
of the process with a real gluon radiated). Note that the integral is not performed 
over the entire phase space, but a lower limit E cut on the gluon energy e g is 
imposed. The reasons for this cut are twofold: first, experimental: only gluons 
with energies greater than a certain threshold can be observed; second, theoretical: 
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when the gluon energy goes to zero, the amplitude A4 7 ' 9 goes to infinity, and the 
cross section itself diverges. This is the well known issue of infrared singularities 
in the radiation of a massless boson. The way to deal with this divergent behavior 
is to note that it is not possible to observe gluons with arbitrarily low energy, no 
matter how good our detector performance. Therefore, the contribution to the 
cross section coming from infinitesimally small gluon energies should be added to 
the cross section for the process without a gluon in the final state []. Thus the 
infrared contributions coming from soft real gluon radiation cancel out the equal 
but opposite-sign infrared quantities which appear in the evaluation of the virtual 
corrections. 

In consequence, the tree level diagrams (Fig. f2- 1|) , the virtual corrections 
diagrams (Fig. |2-4| ) and the real gluon radiation diagrams (Fig. |2-3|) all contribute 
to the NLO cross section for the process with the b, b, W + , W~ final state : 


a 1 (x [ (|A4°| 2 + 2Re[.M 0 (M" 9 )*]) dO 0 + [ \M rg \ 2 dO 0+g (2.8) 

** ** tg^Ecut 

The cancellation of infrared divergences between the real gluon radiation and the 
virtual corrections part is performed by using the phase space slicing method. 
This method amounts to introducing a technical cut-off e for the gluon energy; 


if the gluon energy is smaller than e, the contribution JOT 9 is evaluated in the 
soft gluon approximation. In this approximation the amplitude factorizes to the 


tree-level amplitude times an eikonal factor: 

M " 3 = M " = (2.9) 

(only gluons radiated from the on-shell b quarks contribute to the infrared singu¬ 


larity, as opposed to the case of the on-shell approximation, when gluons radiated 

2 The solution to the problem of infrared divergences is due to Bloch and Nordsieck |24|]; for 
a pedagogical introduction see J23| ; a complete treatment including methods for evaluation the 


soft gluon integral in Eq. 2.1C can be found in [E5l. 
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from the top quarks also give rise to a singular behavior). The phase space also 
factorizes: dfl 0+g = dfl 0 x dfl g ; therefore 


'£ 0 <e 


| M r9 \ 2 dQ 0+g = f \M°{ 2 dO 0 x (c^l 


d 3 * k 


g <e 4:7T£ g 


- 1 WB, 


( 2 . 10 ) 


The term in parentheses (let’s call it Y bb ), which describes the effects of soft gluon 
radiation, can then be computed independently of the tree level amplitude. Since 
it is infrared divergent, Y bb is usually evaluated with the help of some regularization 
procedure (we use mass regularization, which means assuming that the gluon has 
an infinitesimally small mass /r). The same regularization procedure is used to 
compute the virtual corrections diagrams which are infrared divergent (in our case, 
the decay-decay interference diagram |2-4| d) and corrections to the self-energies of 
the b , b quarks). The final result 


a 1 oc / (\M°\ 2 (1 + Y bb ) + 2Re[.A4 0 (Wr 9 )*]) dO 0 + / \M r9 \ 2 dO 0+g 

** ** €<C.6g<I.Ecut 

( 2 . 11 ) 

is finite and independent of the value of the regularization parameter fi. 

Finally, we shall make some comments on the choice of the e and E cut parame¬ 
ters. The E cut is meant to separate between the experimentally observable gluons 
and those gluons soft enough to be undetectable in the experimental setup under 
use. For the detectors planned for an e + e _ linear collider this would mean a value 
for E cut of order 5 to 10 GeV. [] On the other hand, e is an unphysical technical pa¬ 
rameter needed to solve the infrared singularities problem. The only requirement 
on its value is that it is small enough so that the soft gluon approximation works. 
As we shall argue in section 4.4, this means that e should be much smaller than 
the energy of the gluons which contribute to the interference (about 1.5 GeV). 

We therefore shall use for the technical cut e values of order 0.1 GeV. 

3 More about our choice (as well as other cuts intended to select observable gluons) can be 


found in section 3.2 
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2.3 Top width and gauge invariance 


The implementation of finite width for unstable particles is a delicate problem. If 
we use the zero order propagator for the top quark: 

«So(p) = 


- m 0 


the amplitude for our process will have an unphysical non-integrable singularity 
in the region of the phase space where the top propagator becomes resonant. 
Hence the need to regularize this singularity by taking into account the unstable 
particle’s width. However, the introduction of finite width tends to spoil gauge 
invariance. The reason for this is that our theory is gauge invariant in fixed 
order in perturbation theory. The introduction of the width, performed naturally 
through Dyson resummation FI: 


<S — Sq + So(—iE( + ... — 


( 2 . 12 ) 


fi - rriQ - E(fi) 

amounts to taking into account contributions coming from higher orders in per¬ 
turbation theory. Therefore, gauge invariance problems might arise. 

These issues have been first recognized and dealt with in processes where the 
role of unstable particles is played by gauge bosons. There are a variety of methods 
to restore gauge invariance in these cases; the fermion-loop scheme |2fJ is a favored 
one, although it is difficult to implement for complicated processes. For the W 
pair production process at LEP, a simpler prescription is available |y| , [|19 
In the computation of the tree level DPA amplitude: 

Pw+,Pw~) 


M 


DPA ~ 


(Pw+ ~ M'w + iMwTw)(Pw- ~ Mw + iM w T w ) 


(2.13) 


4 Here mo is the bare mass of the particle, while £(^>) stands for the one-particle irreducible 
self-energy. 
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one evaluates the normalized amplitude A4 at the poles (that is, for on-shell 
W’s). The residue M^) is gauge invariant, and the difference between 

this residue and the exact result is obviously non-doubly resonant, therefore it 
can be ignored within the DPA. Thus, in this approach the consequences of hav¬ 
ing off-shell particles are restricted to the denominators of the W gauge bosons 
propagators. Moreover, since the DPA radiative corrections to this process are 
proportional to the tree level amplitude, gauge invariance at NLO is insured also. 

However, the issue of gauge invariance in the W pair production process ap¬ 
pears even at tree level because of the fact that W is a gauge boson. In other 
words, the gauge invariance we are concerned with in Eq. [2.13| is invariance with 
respect to the W gauge. For the top production process, this issue does not arise; 
therefore we can compute the tree level amplitude A4 o with off-shell momenta, 
and the result will be gauge independent because the top quark propagator does 
not depend on any gauge. 

At next to leading order, though, the top quark can radiate a gluon (either real 
or virtual), and the problem of gauge invariance with respect to the gluon gauge 
arises. However, this is a quite different problem from the one discussed above; 
even in the W pair production case, at NLO gauge invariance with respect to the 
gauge of the radiated photon has to be treated separately. For QCD corrections 
to the top pair production case, we shall address this problem for both the case 
of real gluon radiation and virtual corrections in the appropriate chapters below. 

There is one more issue related to the treatment of the unstable particle width. 
The width extracted from the imaginary part of the particle self-energy in Eq. |2. 12 
is a quantity dependent on the top energy - called the running width. In the W 
pair production process, it has been shown that the use of the running width 
generates problems at higher energies [|l^]. Another option is to use the fixed 
width obtained by evaluating the top decay width from the process t —> bW(g), 
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and using this constant value in the top propagator. In the W case, it has been 
found preferable to use the complex-mass scheme , which replaces with the 
complex mass AT 2 = M^—iMwPw not only in the denominator of the propagator, 
but also in the couplings of the particles to the gauge bosons (since these couplings 
can be treated as dependent on the M w mass). In the top production case, the 
couplings do not depend on the top mass; so in our implementation, the top 
mass is replaced by the complex mass m 2 = rri 2 — im t T t (here I\ is the constant 
next-to-leading order top width) in the top propagator denominators only. 



Chapter 3 

Real gluon radiation 


Future high energy lepton colliders — e + e~ and fi + n~ — can provide relatively 
clean environments in which to study top quark physics. Although top production 
cross sections are likely to be lower at these machines than at hadron colliders, the 
color-singlet initial states give lepton machines some advantages. Furthermore, the 
fact that the laboratory and hard process center-of-mass frames coincide greatly 
simplifies the reconstruction of final states. In addition, many of the top quark’s 
couplings, especially those to the photon and Z° boson, can be easily studied 
there. 

The potential for precision studies of top physics at such colliders requires pre¬ 
cision predictions from the theory, beyond leading order in perturbation theory. 
In particular, QCD corrections must be taken into account. One effect of the 
QCD interaction is the radiation of a real gluon. This process can have a siz¬ 
able impact on the analysis of the top production and decay final state variables. 
Jets from radiated gluons can be indistinguishable from quark jets, complicating 
identification of top quark events from reconstruction of top’s decay products. To 
make matters worse, emission may occur in either the top production or decay 
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processes, so that radiated gluons may or may not themselves be products of the 
decay. Subsequent mass measurements can be degraded, not only from misiden- 
tihcation of jets but also from subtle effects such as jet broadening when gluons 
are emitted near other partons. 

In this section we study the effects of real gluons radiated in top quark pro¬ 
duction and decay at e + e _ colliders |^9|. We consider collision energies well above 
the top pair production threshold, so although for definiteness we will refer to 
electrons in the initial state, our parton-level results apply equally well to 
collisions at the same energy. We allow for the top quarks to be off-shell, keeping 
the full width-dependent top propagator and retaining all spin correlations. Gluon 
radiation for off-shell top has been treated previously in the soft gluon approx¬ 
imation [|27], ^]. Here we give an exact treatment for arbitrary gluon energies. 
We study properties of the radiated gluons, top mass reconstruction, and effects 
of interference between production- and decay-stage gluons that can be sensitive 
to the top quark width. 


3.1 Amplitude evaluation 

Our aim is to evaluate the differential cross section for real gluon emission in top 
quark production and decay: 

e + e - 7 *, Z* -»■ tt(g) -> bW + bW~g . (3.1) 

In top events at e + e _ colliders, there are no gluons radiated from the color-singlet 
initial state. Final-state gluon emission can occur in both the production and de¬ 
cay processes, with gluons radiated from the top or bottom quarks (or antiquarks), 
as shown in Figure |2-3| . 

For purposes related to top mass reconstruction, it is desirable to be able to 
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associate the radiated gluon with either the top production process or the decay 
process. While trivial in the on-shell approach, it is not possible to make this 
differentiation with 100% certainty if the top is allowed to be off-shell. However, 
we can use a sensible, ad hoc , definition: if the top is closer to its on-shell mass after 
the gluon has been emitted, call that production stage radiation; otherwise, call 
it decay stage radiation. Thus, emission from the top quark contributes to both 
production- and decay-stage radiation, depending on when the top quark goes on- 
shell. Emission from the b quarks contributes to decay-stage radiation only. The 
separation of these contributions at the amplitude level will be discussed below. 

We compute the exact quantum mechanical amplitudes (also called matrix 
elements) for the diagrams shown in Figure |2-3| with all spin correlations and 


the bottom mass included, using the helicity methods of Kleiss and Stirling [00 


Working at the matrix element, rather than the matrix element squared, level has 
the usual advantages of numerical efficiency, and in our case has the additional 
advantage that we can identify individual contributions as well as interference 
between them. The explicit expressions for the matrix elements are complicated 
and not particularly illuminating, so we do not reproduce them hercQ. A general 
overview including some specific details concerning this computation can be found 
in the appendix. 

We do not assume the top quark to be on-shell; therefore we keep the finite 
top width T f in the top quark propagator and include all interferences between 
diagrams. We use exact kinematics in all parts of the calculation. We do not 
include radiation from hadronic W decays; the W bosons are assumed to decay 
leptonically and we integrate over the decay products in the results presented 
here. In practice, radiative hadronic W decays should probably be taken care of 


1 A FORTRAN program containing the matrix elements can be obtained from the author. 






CHAPTER 3. REAL GLUON RADIATION 


23 


at the level of jet fragmentation simulations; but this requires a separate study. 

3.1.1 Production—decay decomposition 

As mentioned above, calculating at the amplitude level allows us to identify con¬ 
tributions from individual processes and their interferences. We are particularly 
interested in distinguishing between contributions from gluons radiated in the top 
quark production and decay stages. This is directly related to reconstruction of 
the top quark momentum from its decay products, which in an experiment allows 
us both to identify top events and to measure rn t . The presence of gluon radiation 
complicates the reconstruction because the emitted gluon may or may not be part 
of the top decay. If the gluon is not part of the decay, then it is represents a 
correction to top production and should not be included in the top momentum 
reconstruction: 

~ Pt = (Pb + Pw f = Pbw • ( 3 - 2 ) 

(The use of the rs sign in the above equation signifies that the sum squared of the 
4-momenta of the top decay products - for which we will use the shorthand pwb 
- is close to the top mass squared.) If on the other hand the gluon is part of the 
decay, then it should be included in top reconstruction: 

m t ~ Ptg = ( Pb +PW+ Pg ) 2 = P 2 bWg ■ ( 3 ‘ 3 ) 

Being able to make this distinction turns out to be useful for purposes of efficient 
phase-space integration as well. 

Although this production-decay distinction cannot be made absolutely in an 
experiment,":, the various contributions can be separated in the calculation. For 
radiation from the b and b quarks, the assignment is easy: these contributions, 

2 If the interference between processes is large this distinction is not even meaningful. 
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corresponding to the two right-hand diagrams in Fig. |2-3| , are clearly part of the 
top quark decay. However, as noted above, gluon emission from the top quark (or 
antiquark) contributes to both the production and decay stages; which is which 
depends on whether the top was closer to its mass shell before or after emitting 
the gluon. This condition corresponds to which of the two propagators from the 
top that emitted the gluon is numerically larger. 

We can make the separation in our calculation as follows . For definiteness, 
we consider gluon emission from the top quark, shown in the upper left diagram 
in Fig. |2-3| . The matrix element for this diagram contains propagators for the top 
quark both before and after it radiates the gluon. The matrix element therefore 
contains the factors 


M t oc 


\Pwbg - m t + irrit^t J \Pwb ~ m t + t, 
The right-hand side can be rearranged to give 


(3.4) 


1 1 

Pw h ~ rrif + irrifTt pf Vbg - rrif + irn t V t 

This separates the production and decay contributions to the matrix element. 
The first term in parentheses contains a propagator that peaks when p^ Vb = mf, 
which corresponds to the condition for production stage, as in Eq. |3.2| . Similarly, 
the second term peaks for pf Vbg = rrif , which corresponds to decay emission as in 

Eq. |[3|. 

The complete amplitude in Eq. [2.4| can now be rewritten schematically as Q 



M t oc 


2 pwb * P a 


AA. ' AAprod T AAtdecay T AAtdecay • 


(3.6) 


The cross section, obtained from taking the absolute square of M rg , then con¬ 
tains separate production and decay contributions, from |-M pr . 0 d| 2 and \Mtdecay\ 2 , 


Expressions for the partial amplitudes can be found in Appendix B. 
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\Aiuecay\ 2 i respectively. It also contains cross terms representing the interferences, 
which in principle confound the separation but in practice are quite small. 

The interference terms are interesting in their own right, although not for top 
reconstruction. In particular, the interference between production- and decay- 
stage radiation can be sensitive to the top quark width T t |27|, |28|, which is 1.42 
GeV in the Standard Model at 0(a s ) 0 . The interference between for example 
the two propagators shown in Eq. |3l] can be thought of as giving rise to two 
overlapping Breit-Wigner resonances. The peaks are separated roughly by the 
gluon energy, and each curve has width T t . Therefore when the gluon energy 
becomes comparable to the top width, the two Breit-Wigners overlap and there 
can be substantial interference. In contrast, if the gluon energy is much larger 
than T t , the overlap — and hence the interference — is negligible. Therefore the 
amount of interference serves as a measure of the top width. We will explore this 
more below. 

3.1.2 Gauge invariance 

Since gauge bosons contribute to the processes analyzed in this thesis, it is natural 
to consider the issue of gauge invariance. In order to perform our computations, it 
is necessary to choose some gauge for the bosons involved (for example, the propa¬ 
gators for the 7 , Zq bosons are computed the Feynman gauge, while the definition 
of the gluon polarization vector depends on a gauge parameter as described in 
the appendix). While our choice of gauge may vary, the final results should be 
independent of this choice. 

Due to the fact that the electrons in the initial state are considered massless, 
there is no amplitude dependence on the gauge in which we evaluate the 7 , Z 0 
propagators. However, this is not true anymore for the case of the gluon gauge. 
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This has to do with the fact that the top quarks can be off-shell; it is an easy 
exercise showing that amplitudes in the on-shell approximation are independent 
of the gluon gauge. More precisely, the reason that gauge invariance is lost is that 
the physical process we actually compute is: 


bW + b W~g 


(3.7) 


and to obtain a gauge invariant result we should take into account all the Feynman 
diagrams contributing to this final state. Since our computation includes only 
diagrams with two intermediate top quarks, the final result is not strictly gauge 
invariant. 

ffowever, the diagrams we take into consideration are the only ones from the set 
contributing to (|3.7|) which have a doubly resonant structure. So, if one considers 
the gauge invariant result as being a sum of doubly-resonant, singly resonant 
and nonresonant terms, the only diagrams which can contribute doubly resonant 
terms are the diagrams containing two top quarks in Figure |2-3| . Therefore, the 
amplitude in ( |3.6|) differs from the gauge invariant result by singly or non-resonant 
terms only. This means that by subtracting non-doubly resonant terms from A I tot 
(which is allowed in DPA) we can obtain a gauge invariant answer. 

One way to perform these subtractions is as follows. Consider the diagram 
where the gluon is radiated by the top quark (the first diagram in Fig. |2-3| ). The 
contribution of this diagram to the production amplitude is: 



By commuting /k g to the right, this can be written 


M (t) _ _J_ ( j>wb + mt)(2e g • p W b+ 4 g fc) - ( p 2 Wb - m t 2 ) 4 g 



...v(b) . 
(3.9) 
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The term which breaks gauge invariance here is the one proportional to (pl Vb — 
m 2 ). This is a non-resonant term, regardless of the gluon being radiated in the 
production or decay stage (in other words, regardless of pyy b ~ mf or Pyy bg ~ mf); 
therefore, in keeping with the approximation used, we can neglect it. 

A similar analysis works for the contribution of this diagram to the top decay 
amplitude. Starting with the expression: 


M 


(t) 

tdecay 


u{b) 4w ^ + 4g (ftwbg + m t )... v(b) (3.10) 


2 kpt ' ' ■ Pwbg - m t 

by commuting / g to the left, we obtain: 


M 


(t) 

tdecay 


2kp, 


■u(b ) Av 


MPwbg - rn 2 t ) - (2e 9 • p Wbg - Jk tg)(ftwb g + m t ) 


PlVbg ~ m t 


... v(b) 

(3.11) 

and in this case we drop the term proportional to (py/bg ~ m t) ■ Finally, the 
amplitudes corresponding to the diagram in which the gluon originates from the 
t can be computed in the same manner. The final result is gauge invariant, and 
differs from the exact result by non-doubly-resonant terms, as we have shown. 

We have implemented the above computation in the Monte Carlo program, 
and have checked numerically that the difference between the gauge invariant 
result and the exact result is very small (of order 0.01% of the total cross section, 
and order 1% with respect to the interference terms). This indicates that the 
other non-doubly resonant contributions (coming from diagrams with a single top 
or none) are also small; a more detailed study is in progress. 

Finally, we note that this method for restoring gauge invariance is not unique. 
We could, for example, have instead replaced the top mass in the top propagator 
numerator with the invariant masses: \Jpwb i n the production amplitude, and 
\JPwbg i n the decay amplitudes. The result obtained with this method is also 
gauge invariant, and also differs from the exact result by non-resonant terms. 
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3.1.3 Monte Carlo and phase space integration 

The integration over the final state phase space to obtain the cross section involves 
an integrand that contains multiple Breit-Wigner peaks from the top quark prop¬ 
agators as well as infrared singularities when the gluon energy becomes small. 
Even with cuts on E g , the rapid variation of the integrand can spoil the integra¬ 
tion procedure. To eliminate this problem, we tailor the momentum generator to 
the production of a gluon in association with two massive particles ( 7 *, Z* —»ttg 
or t —> bWg ). The multiple Breit-Wigner peaks are taken into account by using 
a multi-channel approach that integrates separately over the individual produc¬ 
tion and decay contributions; the Breit-Wigner behavior is smoothed out in the 
phase space generation. The interference terms, which have products of Breit- 
Wigners that peak in different places, much like in Eq. |3.4| , are integrated using a 
combination of the three main channels. 


3.2 Numerical Results 

In this section we show results of the numerical calculation described above. We 
present the cross section for bbW + W~g production in e + e _ collisions at a 500 
GeV center-of-mass energy, with a few exceptions which are clearly identified. 
The calculation is entirely at the parton level, and we do not include initial state 
radiation, beam energy spread, or beamstrahlung. We use the following numerical 
values of parameters: m t = 175 GeV, mb = 5 GeV, Mw = 80 GeV, T t = 
1.42 GeV, and a s = 0.1. Note that for the results presented in this section a s 
appears simply as an overall factor, because all of our events contain a gluon. 

Unless otherwise indicated, we use the following cuts. We require E g > 5 GeV 
to eliminate the infrared singularity and because we intend for the gluon to be 
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detectable. In addition we wish the gluon to be separable from the b and b 
quarks; this is implemented via the requirement rrib g , rrii g > 10 GeV, which we 
shall identify below as u rrib g cuts.” (Separation could also be achieved with a 
cut on the gluon’s transverse energy Et with respect to the b or b\ the choice 
makes little difference in the resulting distributions.) In order to make sure that 
we do not get contributions to our results from regions of the phase space where 
non-doubly-resonant diagrams might be important, we require 

160 GeV < rribw < 190 GeV or 160 GeV < m b w g < 190 GeV (3.12) 

and the same thing for the b. These conditions will be identified as “m t cuts” 
below. 

3.2.1 Characteristics of the gluon radiation 

We begin with the relative contributions of production- and decay-stage radiation 
to the total cross section. Figure |3-1| shows the fraction of the total cross section 
due to production stage emission in events with an extra gluon, as a function 
of the minimum energy of the gluon. This figure contains no cuts besides that 
for gluon energy and is simply meant to illustrate how radiation is apportioned 
in top production and decay for different center-of-nrass energies; the solid line 
corresponds to c.m. energy 1 TeV, and the dashed line is for 500 GeV. Both curves 
fall off as the minimum gluon energy increases; this reflects the decrease in phase 
space for gluons radiated in the production stage. We see that the production 
fraction is always higher at 1 TeV collision energy than at 500 GeV. This too 
reflects phase space — for a given gluon energy there is more phase space available 
to produce gluons in association with top pairs at the higher c.m. energy. However 
both fractions remain below 0.5; decay-stage radiation always dominates at these 
energies. 
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Figure 3-1: The fraction of gluon emissions radiated in the production stage, as a 
function of minimum gluon energy, for center-of-mass energy 1 TeV (solid line) and 500 
GeV (dashed line), with no cuts. 

Figure |3-2| shows for a 500 GeV center-of-mass energy the effect on the pro¬ 
duction fraction of separation cuts between the gluon and b quarks. The dashed 
line shows the fraction with no cuts. The dotted line corresponds to requiring 
that the transverse energy of the gluon with respect to the b and b — which we 
denote E T (g , b) — be greater than 3 GeV. The solid line corresponds to the cut 
rribg > 10 GeV where the b can either be a quark or antiquark. The effect of both 
of these cuts is to eliminate gluons that are soft and/or close to one of the bottom 
quarks; since these contributions tend to come from decay-stage radiation, their 
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Figure 3-2: The fraction of gluon emissions radiated in the production stage, as a 
function of minimum gluon energy, for center-of-mass energy 500 GeV, with no cuts 
(dashed line), Er(g,b ) > 3 GeV (dotted line), and m&g > 10 GeV (solid line). 

effect is to increase the fraction of production-stage radiation. If the b were mass¬ 
less there would be a collinear singularity in the decay contribution; this does not 
happen in our case but the decay distribution still peaks when the 6-quark-gluon 
angle is small. The effects of both cuts become smaller with increasing gluon 
energy. 

Figure |3-3| shows the total gluon energy spectrum for a collision energy of 
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Figure 3-3: The spectrum of radiated gluons as a function of gluon energy in GeV for 
center-of-mass energy 500 GeV, with m^g and cuts (see text). Dashed histogram: 
production-stage radiation. Dotted histogram: decay-stage radiation. Solid histogram: 
total. 

500 GeV along with its decomposition into production (dashed histogram) and 
decay (dotted histogram) contributions. The interferences between the two are 
negligible and are not shown; this will be true for all subsequent figures until we 
consider the interference explicitly. Included in this figure are the m t and m^g 
cuts discussed above. As indicated in the previous figures, radiation from the top 
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Figure 3-4: The top invariant mass spectrum without (left) and with (right) the gluon 
momentum included, for center-of-mass energy 500 GeV, with m bg cuts and E g > 5GeV. 

decays dominates. Otherwise the spectra are not vastly different; both exhibit the 
rise at low energies due to the infrared singularity characteristic of gluon emission, 
and both fall off at high energies as phase space runs out. 

3.2.2 Mass Reconstruction 

We now turn to the question of top mass reconstruction in the case when there 
is a gluon in the final state (more about measuring the top mass can be found 
in the next chapter). Figure |3-4| shows top invariant mass distributions with and 
without the extra gluon included; the first plot shows the distribution in rn b w and 
the second shows m bWg . We have imposed rn bg cuts and required E g > 5 GeV. In 
both cases there is a clear peak at the correct value of rn t . Note that the peak in 
the first plot contains the production contribution as expected, but the radiative 
decay part contributes as well. This is because even for decay-stage radiation, only 
one of the produced t quarks decays radiatively; the other still has p\ = pf )W and 
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Figure 3-5: The distribution in the angle between the gluon and the b quark for 
center-of-mass energy 500 GeV, with E g > 5 GeV. The various contributions are as 
described in the text. 

therefore contributes to the m^w peak. The long tails in the two distributions are 
from misassignments of the gluons. In the left-hand plot, where the gluon is not 
included in the reconstruction, we see a low-side tail due to events where the gluon 
was radiated in the decay but was not included in the reconstruction. Similarly, 
in the right-hand plot we see a high-side tail cine to events where the gluon was 
radiated in association with production, and was included when it should not have 
been. 
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We wish to define a single distribution for the top mass that combines both 
types of events yet omits wrong combinations as much as possible. One possibility 
is to cut on the angle between the gluon and the b quark, whose distribution we 
show in Figure f >-5| . This is motivated by the fact that gluons radiated from the 
b quarks are always part of the decay, and such gluons tend to be emitted close 
to the b quark direction. As we have mentioned, the mass of the b quark prevents 
a collinear singularity, but the gluon distribution still peaks close to the b , as 
can be seen in the figure. Because we wish to define cuts that give a narrow 
top invariant mass distribution, the distribution in 6 bg is decomposed into various 
invariant mass regions. (Here we refer to the b quark only, and not the b.) Using 
the variables rh t = m bW +, fh tg = m bW + g , fht = m lw- an d = m lw~g we define 
four types of events: 

type 1 : 172 GeV < fh tg ,fht < 178 GeV (vertical hatching) 
type 2 : 172 GeV < m t , ffiig < 178 GeV (horizontal hatching) 
type 3 : 172 GeV < m t ,fhi < 178 GeV (cross hatching) 
type 4 : any other event (no hatching) 

Type 1 events are dominated by contributions from radiative t decays, and we can 
see that they do tend towards the b quark direction. Type 2 events (horizontal 
hatching) are in turn dominated by radiative t decays; gluons in this case tend 
to cluster near the b direction, and since the b and b tend to appear in opposite 
hemispheres, type 2 gluons are mostly found at large angles to the b. Events of type 
3 (cross hatching) are mostly pro duct ion-stage contributions; their distribution is 
more or less uniform, independent of the b quark direction. Finally, events of 
type 4 (no hatching) get contributions from both production and decay, with no 
compelling evidence for one over the other. 

Using this figure we can make the following conventions: 
if Q bg <0.7 rad, assign gluon to t decay 
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Figure 3-6: The top invariant mass spectrum with 6-gluon angle selection criteria 
(dotted histogram), for center-of-mass energy 500 GeV, minimum gluon energy 5 GeV, 
and rribg cuts. The solid curve and histogram show the effects of energy smearing. 

if 6} ig <0.7 rad, assign gluon to t decay 

if dbg, dig > 1 rad, assign gluon to tt production. 

With these cuts on the proximity of the gluon to the b quark, we construct the 
top mass distribution presented in the dotted histogram in Figure fi-6| . 

Of course an important reason the cuts are so effective is that we work at the 
parton level. The experimentalists do not have that luxury, and, as one would 
expect, hadronization and detector effects are likely to cloud the picture. The solid 











CHAPTER 3. REAL GLUON RADIATION 


37 


histogram in Fig. |3-6| shows the mass distribution after including energy smearing; 
the solid curve is a Breit-Wigner fit. The spread in the measured momenta of the 
final state particles is parametrized by Gaussians with widths o = 0 Ay/E for 
quarks and gluon, and a = 0.15 y/E for the W's. We see that the central value 
does not shift, but the distribution becomes significantly wider. 

These results are meant to give an indication of the effects of hard gluon radi¬ 
ation on mass reconstruction and how they might be dealt with. Other variables 
to consider in choosing the cuts are m \, g , the transverse energy of the gluon with 
respect to the b or b, or some combination of energies and angles as defined in the 
various algorithms used in jet definitions for e + e _ colliders. At tree level and with 
partons only, the exact choice is not very important. We will revisit the question 
in more detail when we include virtual corrections in a full NLO calculation. 


3.2.3 Interference and Sensitivity to F t 

Finally, we return to the subject of interference. As mentioned above, the inter¬ 
ference between the production- and decay-stage radiation can be substantial for 
gluon energies close to the total width of the top quark the interference is there¬ 
fore sensitive to the value of T t . However, because this interference is in general 
small, we need to find regions of phase space where it is enhanced. This question 
was considered in |28| in the soft gluon approximation^, where it was found that 
the interference was enhanced when there was a large angular separation between 
the t quarks and their daughter 6’s. 


Here we examine whether the result of |28[] , which considered a fixed final-state 
configuration, survives the exact calculation and phase space integration. Figure 
|3-7| shows that it does. There we plot the distribution in the angle between the 


4 See also [27, 3^] 
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0 tg (degrees) 


Figure 3-7: The distribution in angle between the top quark and the gluon for gluon 
energies from 5 to 10 GeV, cos Otb, cos 6^ < 0.9, mt cuts, and 750 GeV collision en¬ 
ergy. The upper solid histogram is the total and the other histograms represent the 
individual contributions: dotted: decay; dashed: production; dot-dashed: decay-decay 
interference; solid: production-decay interference. 


emitted gluon and the top quark for gluon energies between 5 and 10 GeV and with 
cos 6 t b < 0.9 and m t cuts. The center-of-mass energy is 750 GeV. This c.m. energy 
is chosen because for there to be significant interference between production and 
decay-stage radiation, both contributions must be sizable. At 500 GeV, we see 


from Figs. 3-1 and 3-2 that the production contribution is suppressed compared to 
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that from decay; as a result, the interference is very small. Increasing the energy 
increases the pro duct ion-stage contribution. We note that the distributions at 750 
GeV and 1 TeV do not differ substantially. 

The histograms in Fig. |3-7| show the decomposition into the various contri¬ 
butions. The production-stage radiation is shown as a dashed histogram; we see 
that it reaches its largest values at relatively small and large angles. Small angles 
correspond to the t direction, and large angles more or less to the t direction, since 
for the small gluon energies of interest here, the t and t are nearly back-to-back. 
The dotted histogram represents the decay-stage contribution; it dominates the 
cross section and peaks in the same region as the production contribution. This 
leads to substantial production-decay interference, shown as the negative solid 
histogram. This interference is destructive, so that it serves to suppress the total 
cross section, shown as the positive solid histogram. This effect would be en¬ 
hanced if we lowered the gluon energies to values closer to T t , but jets from very 
low energy gluons are not likely to be observable, so we cut off the gluon energy 
at 5 GeV. Finally, interference between the emissions in the t and t decays are 
shown as a dot-dashed histogram, but as there is very little overlap between the 
two phase space regions even with these cuts, this contribution is negligible. 

The cuts we have used are fairly generic; we can further enhance the interfer¬ 
ence terms with a judicious choice of additional cuts. If we examine their behavior 
in more detail in various regions of phase space, we find that the sign of the in¬ 
terference terms depends on the value of the invariant mass of the top quark. 
Since we integrate over this mass, we get cancellations (a similar effect ensures 
cancellations of non-factorizable corrections in inclusive quantities). 

Consider the interference between radiation in the production stage and the 
top decay stage. The product of the two Breit-Wigner peaks is proportional to 
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Figure 3-8: As in Fig. |3-7|, with the addition of the cuts given in Eqs. |3. 14.3.15. 


the factor 

ft = (Pwb ~ m 2 t ) * ( p 2 Wbg - m 2 t ) + m\T 2 (3.13) 

This factor will multiply a quantity which, upon integration over angles, is nega¬ 
tive. Therefore, for invariant mass values such that f t is positive, the interference 
terms are negative, while for negative ft, the interference terms are positive. We 
can impose cuts that take advantage of this: if we require the invariant masses to 
satisfy 


ft > 0 , if % < 


(3.14) 
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Figure 3-9: The distribution in angle between the top quark and the gluon for gluon 
energies from 5 to 10 GeV, cos Otb, cos 6^ < 0.9, and 750 GeV collision energy. The 
histograms correspond to different values of the top width T^: dot-dashed: 0.1 GeV; 
solid: 1.42 GeV (SM); dashed: 5. GeV; dotted: 20 GeV. 

ft > o, if e bg > e h , (3.i5) 

we obtain the distribution shown in Figure [3-8[ . The interference effects are en¬ 
hanced, though at the cost of a substantial decrease in cross section. 

Because the production-decay interference is destructive, increasing the top 
width would further suppress the total distribution. The height of the peaks, 
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then, is sensitive to the value of IV This is illustrated in Figure |3-9| , which shows 
the cross section (without the cuts of Eqs. j3.14| , |3.15D for different values of the 
top width]], ranging from 0.1 GeV to 20 GeV. The SM case (T t = 1.42 GeV) is 
shown as a solid line. It is interesting to note that in the context of perturbative 
gluon radiation, the SM top width is actually a small quantity. There are sev¬ 
eral other points to note. In principle, this sensitivity to T t gives us a method 
to measure the top quark’s total width, independent of decay mode, above the 
top production threshold. Although in practice statistics would surely limit the 
possible precision of such a measurement, the total top width is not so easy to 
measure directly by any means. Furthermore, the effects illustrated here arise 
from simple quantum-mechanical interference, and finding experimental evidence 
for interference between the radiation at the various stages is an interesting goal 
by itself. 


5 The histograms here are scaled so that they all would have the same normalization in the 
absence of interference effects. Without this rescaling, changing the width while keeping the 
tbW coupling fixed changes the total cross section, which behaves like 1/T^ for small r t . 








Chapter 4 

Virtual corrections 


The subject of this chapter is the computation of virtual corrections to the process: 

e+e" -► 1t —> b W + b W~ (4.1) 

Although real gluon radiation plays an important role in top production and de¬ 
cay, the four particle final state in Eq. E3 will be the predominant signal observed 
experimentally. Therefore, it is of utmost interest to have good theoretical pre¬ 
dictions for this process, and this means that the virtual corrections to the tree 
level amplitude should be computed. 

While in studies of top quark production near threshold most of the informa¬ 
tion is extracted from the shape of the total cross section as a function of the 
beam energy, above threshold the information will be extracted from the analysis 
of kinematical variables of the decay products. Consequently, an analytic compu¬ 
tation is much less useful in this case (and much less feasible, too), so the approach 
we will use to perform our computations will be Monte Carlo simulation. 

One of the most important quantities to measure at future collider experi¬ 
ments is the top quark mass. A precise measurement of this parameter will allow 

stringent tests of the Standard Model of particle physics (or give information on 
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the supersymmetric model, if SUSY has been discovered by that time). Indeed, 
it is well known that, at the present time, the measured top quark mass together 
with values for the W boson mass and other electroweak parameters can be used 
to constrain the Higgs mass in the Standard Model (these constraints are not very 
strong so far, but they will be greatly enhanced by improvements in the top quark 
mass measurements p2|). 

One way to measure the top mass at energies above the production threshold 
is by determining the position of the peak in the distribution of the invariant mass 
of the top decay products ^( pw + Pb ) 1 2 (this is the way the top mass has been 
measured at the Tevatron). The quantity thus obtained is called the pole mass; 
the threshold studies will measure a quantity called the threshold mass (closely 
related to the MS mass) []. A theoretical relation exists among these quantities; 
hence a measurement of both of them will allow a test of our understanding of 
the underlying theory. 

The extraction of the top quark mass from the bW invariant mass distribution 
is the most natural and straightforward method; therefore, there are a lot of studies 
in this area (and we shall present some results of relevance to this case, too). 
However, it suffers from a fundamental limitation on the precision it can achieve: 
since the physical observable is not a color singlet, there will be uncertainties 
of order Aqcd in the measured quantities (what this actually means is that we 
cannot measure the momenta of the b jet - of any jet generated by a parton which 
has color, for that matter- with a precision greater than Aqcd)- Therefore, the 
minimum error achievable through this method is about 200 MeV. In order to 
circumvent this limit, it might be desirable to try to extract the top mass from 

1 Since a quark is not a stable particle, its mass is an elusive quantity; there are several 

theoretical definitions for the mass of the top quark; among these, the most common are the 
pole mass and the MS mass. 
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distributions which involve color singlet physical observables; for example, end 
points in distributions of quantities like Pw(Pb+Pb)i or the line shape of the decay 
leptons energies However, since these methods are quite involved, and would 
also require a good understanding of hadronization and detector issues, we will 
not discuss them further. 

This kind of precision measurement obviously cannot be performed without 
a good theoretical understanding of the underlying process. Hence the need to 
compute the QCD corrections to the top production and decay process. A quick 
analysis shows that these corrections are indeed important, ft is well known that 
the total NLO QCD corrections to the production of a pair of massless quarks is 
ccs/ttx the Born cross section. For the top production cross section, this would 
correspond to the high energy limit (beam energies of order TeV). On the other 
hand, near the top production threshold, the QCD corrections are quite large; 
from the results presented in JT0|. we have: 


as j3 —> 0 

where (3 is the top quark speed []. This indicates that the NLO corrections to top 
production are rather large for CM energies smaller that 1 TeV; as we shall see 
below, they are about 20% of the Born cross section at 500 GeV, and 7% at 1 
TeV. 

In the following, we shall present a NLO computation for the top production 
and decay process which takes into account not only corrections to the produc¬ 
tion and decay subprocesses, but also interference contributions. The framework is 
2 Near the production threshold, fixed order NLO computation does not work anymore; rather, 
we need to perform a resummation of the large log/3 logarithms. Following the suggestion 
advanced in jn]], we can assume that the NLO computation gives reasonable results if a s /f3 < 
1/4, which means that the collision center-of-mass energy should be bigger than about 380 GeV. 


ale e 
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Born 
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that of the double pole approximation, as described in Chapter 2. The contribut¬ 
ing amplitudes (Eq. [2 . 5| ) are computed using off-shell momenta, and numerical 
integrations for the cross section are performed over the full off-shell phase space. 
This approach allows for further development of our Monte Carlo (taking into ac¬ 
count singly resonant terms, for example) to be performed in an easy and natural 
way. We shall also discuss an alternative computational method (on-shell DPA) 
and compare the results with those obtained in our approach. Some results on the 
total cross section for top production and decay and the impact of interference 
corrections on top mass reconstruction are presented at the end. 


4.1 Resonant structure of partial amplitudes 


In this section, we shall discuss the evaluation of NLO amplitudes contributing 


to the process |4. 1| . Some of the relevant Feynman diagrams are presented in 
Figure 2-4. As mentioned before, these diagrams can be split into two classes: 
(i) vertex and fermion self-energy diagrams, and (ii) interference diagrams. In 
the following, we shall look at each of these types separately, starting with the 
interference diagrams. 

However, let’s first consider some general properties of these amplitudes. In 
keeping with the approximation used (DPA), all of them have a doubly resonant 
behavior in the phase space region where the invariant mass of the Wb pair is close 
to the top mass. However, the exact type of behavior is different from diagram 
to diagram; while for the zero order amplitude the resonant behavior is of type 
pole x pole: 

. . 1 1 
AA 0 




j - m; 


pi - m P i 

for the interference amplitudes the resonant behavior is formally of type pole x 
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logarithm ; for example: 

Mbt ~ lo gfe 2 - vA— 2 

Pt ~ 

The replacement of a po/e-type singularity with a Zog-type singularity can be 
traced back to the integration over the momentum of the gluon in the virtual 
loop. Moreover, this kind of term does not appear only in the purely interference 
diagrams, but also in the vertex and self-energy diagrams. The amplitude corre¬ 
sponding to the top decay vertex correction diagram, for example, can be written 
as a sum of terms: 


Mtb = 


2 — o 2 — 2 

pt - Pi ~ m i 


Mtdec + log (Pt - m 2 t ) 


l2 M' tb 


Pi - m i 

where the terms with the double pole structure can be thought of as being the 
correction to the top decay process, and the term with a pole x log structure 
contributes to interference between the production and top decay process. 

Another property of the amplitudes for the interference diagrams is that, in 
the DPA, they are proportional to the tree level amplitude A4q. We say that, in 
the DPA, the interference amplitudes factorize. This is also true for the case of 
the W pair production process (|^j, |H]); moreover, for this process it has been 
shown that all the interference contributions (including those coming from vertex 
correction diagrams, for example) factorize. Therefore, the total interference can 
be written as a scalar factor times the Born amplitude, which makes the evaluation 
of these contributions really easy. However, this is not true in our case; for the top 
production and decay process, there are interference contributions coming from 
the vertex and top self-energy correction diagrams which are not proportional to 
the tree level amplitude, as we shall see below. 
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4.1.1 Interference diagrams 

In this section we will examine the behavior of interference diagrams. Consider, 
for example, the t — b interference diagram. (Fig. 2c). The amplitude associated 
with this diagram is: 




(4.2) 


where the Feynman gauge is used for the gluon propagator. 

The evaluation of this amplitude is obviously quite a difficult task. However, 
the only terms of interest to us in DPA are those which have resonances at the top 
and antitop quark propagator poles. This simplifies our task greatly. The doubly 
resonant terms can be extracted with the help of the following observation: if the 
virtual gluon in the loop is hard, then the quantity in brackets does not have any 
singularity, and the overall resonant structure for this diagram is given only by 
the pole due to the antitop propagator: .M&j(hard gluon) oc l/(pf — fhf). This 
means that any doubly resonant terms contribution to A4 b t are entirely due to 
soft virtual gluons. Therefore, we can neglect the /k terms in the numerator of 
(|4.2|). Following |I9[, we shall call this approximation the extended soft gluon 
approximation (ESGA) []. 

With the help of the transformations: 


Y (A + m b ) = (- A + m b ) + 2 Pb -*■ 2 Pb 


(4.3) 


(- A + m) = 7 m (A + m) - Zptu, -*■ -2 Pin 


3 In the standard soft gluon approximation, fc 2 terms in the denominator of top quark propa¬ 


gators would also be neglected; we do not do this here for computational reasons (see also [ p8| , 
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(the term + on the second line is neglected, since it would lead to a singly 
resonant contribution), the following result is obtained for the amplitude ( |4.2|) : 


ry 

M b t(DPA + ESC A) =-± Mo * (-4 p b pi)(p 2 t - m?) * 

47T 

[ — — - - - - - -_- - -- (4.4) 

J iir 2 k 2 +ie k 2 — 2kp b (pt — k) 2 — m\ {pi + k) 2 — m 2 

The result is proportional to the tree level amplitude - in the DPA, the virtual 

corrections due to interference factorize. The proportionality factor includes the 

scalar four point function (the integral on the second line of Eq. |4.4|) T>^ = 

'D°{~Pbi^PtiPt-> 0, rn b , fh t , fh t ). 1] 


What can we tell about the singular behavior of the DPA amplitude in ( [4.4] )? 
Apparently, the result for M b i has a single pole at p\ = m 2 (the other pole being 
canceled by the multiplicative term p 2 — mf). However, if the top (or antitop) goes 
on-shell, the T>° function acquires an infrared singularity (in the zero top width 
limit; this singularity is regularized by the top width). Since the infrared singular 
type terms have a logarithmic structure (this can also be reasoned from power 
counting arguments), this indicates that D° has the following behavior close to 
the top resonances: 


V° bl ~ cqlog(^ - m]) + a 2 log(^ - m\) 

4 For the scalar one-loop integrals appearing here we use the following notation: 


T> 0 (pi,p2,P3,mo,m.i,rn2,rn 3 ) = 

£°(pi,.P2,P3,P4,too, wi,rn 2 ,m 3 ,TO 4 ) 
with the denominators : 


d A k 1 
in 2 N 0 W N 2 N 3 

r d 4 k 1 

in 2 N 0 Ni N 2 N 3 N a 


(4.5) 


N 0 = k 2 - ml + it, Ni = (fc + pi ) 2 - m 2 + it, 


i= 
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Here cq and cq are terms which are finite when either the top or antitop quark go 
on-shell. 

Formally, then, the overall resonant behavior of the interference amplitude 
A4bi in DPA is of type pole x logarithm : 

M b i ~ Mo log (p 2 - rh 2 t ) 1 

Pi ~ mf 

rather than pole x pole, as it is for the corrections to production or decay subpro¬ 
cesses. 

Using the same techniques, similar results are easily obtained for the other two 
interference diagrams. In the soft gluon approximation (and DPA): 

M t i(DPA + ESGA) = Az Mo * (-4 m )(p| - m 2 ) P? 5 (4.6) 

M bb {DPA + ESGA) = -A Mo * (-4 p b p- b )(p 2 - m 2 )(p 2 - m 2 ) £& (4.7) 

where V { I = V°(-p- b , -pi, p tl 0, m b , rh t , rh t ) and 


£ bb = £°(-p bl -pt,Pi,p b ,g,m b ,mt,m t ,m b ) (4.8) 

is the scalar five point function (here p is the infinitesimally small gluon mass 
needed for the regularization of infrared divergent behavior of £ bb ). 

We end this section with some comments on the numerical magnitude of in¬ 
terference terms. Since the resonant behavior of these terms is of pole x log type, 
it might be expected that they are less important numerically that the double 
pole terms. However, analytic expressions for the V 0 function ([[L9|, |2lH ) 

show that, although the coefficients Gq,a 2 in |4.5| are finite when one of the top 
or antitop quark goes on shell, they will diverge when both particles go on-shell 
simultaneously: 


1 

cii (Pt ~ rht.) + C2i ( P 2 ~ rh 2 t ) 

Therefore, the leading logarithms in the scalar 4 and 5-point functions will be 
enhanced by factors of order mt/T t near the top, antitop quark mass resonances. 
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Figure 4-1: General vertex correction diagram. 


4.1.2 Vertex corrections 


The results for the interference diagrams are completely analogous to results ob¬ 
tained in the W pair production computation. However, for the off-shell vertex 
and self-energy corrections diagrams, the results in the top case are different. 
Consider for example, the amplitude for the general vertex correction in Figure 
041 




47T 


d 4 k 1 

77T 2 k 2 


7 


fa- k + 

(pi ~ k ) 2 — fh\ 


r(c v + c A 7 5 ) 


- fa- /k + m 2 

(. P2 + k ) 2 - m 2 ^ v 


Upon evaluation (and keeping only the vector part) the result can be written in 
terms of eight form factors, each of them multiplying a different tensor quantity: 


(y 

/CV [faF 2 + (A - m 1 )'fF A + 7 m (- fa ~ m 2 )F e + 

47r 

(fa - mi)fa(- fa - m 2 )F 8 +plF 1 + (fa- mi)p^F 8 + ...] (4.9) 

(expressions for the scalar form factors Fi,... , F 8 can be found in the appendix). 
In the on-shell case, only the F 2 (electric dipole) and F\ (magnetic dipole mo¬ 
mentum) form factors contribute, ft might be expected that in the double pole 
approximation we can drop the other terms, too, since they have a zero at fa = m\ 
(or fa = —m 2 ) which will cancel one pole (or both) in the amplitude. However, 
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the form factors themselves may have a resonant structure when the particles go 
on shell. 

Consider the top decay vertex correction. In this case, p\ —> Pb,P 2 —> ~Pt, and 
only four terms survive in Eq. |4.9| ; the corresponding form factors contain terms 
which are proportional to the scalar three point function: 



i = 1,2,5,6 


which has a logarithmic resonant behavior: 


Ctb ~ log (p 2 t -m t 2 ). 


Therefore, the contribution of i = 5, 6 terms to the top decay vertex correction 
is doubly resonant, although of type pole x log rather than double pole. Because 


these logarithms are not multiplied by large factors (as in the case of the interfer¬ 


ence diagrams), we can expect these terms to be numerically small; for consistency 
reasons it is still desirable to include them in the final result. 

Similar results are obtained for the correction to the antitop decay vertex (we 
keep the i — 1, 2, 3,4 terms in this case). In the case of the t — t vertex, though, 
both fermions are off-shell; as a consequence, there are no resonant logarithms 
when either the top or antitop quark goes on-shell, and we keep only the i — 1,2 
terms. 

It follows that in the general expression ([0|) we have to keep the terms which 
contain F\ to Fq (we can drop the F-j and Fg terms), and we don’t have factor¬ 
ization anymore. This is different from what happens in the W pair production 
process, where in DPA factorization holds even in the off-shell case. This differ¬ 
ence is due to the fact that in our process the intermediate particles are fermions, 
and not bosons. 
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Figure 4-2: Terms contributing to the renormalized vertex; the dots represent 
counterterm insertions. 


4.1.3 Renormalization and fermion self-energy 

Since we are performing a next-to-leading order computation, we have to deal with 
the issue of ultraviolet divergences and renormalization. We use the counterterm 
method (for a presentation of the essential features see for example ||23|). What 
renormalization amounts to in our case is replacing the bare vertex correction in 
Figure |4-1| (which is UV divergent) with a finite renormalized vertex correction: 




STL 


ar e „ = + + -— 

2 p- L - mi 


rM + F rM - T~ - 1 

2 - pi ~ rn-2 


2 ) 

(4.10) 


in which we have included also the contributions of the fermion self-energy dia¬ 


gram. The diagrams corresponding to separate terms in Eq. |4.10| are presented 
in Figure |4-2| . 

The first two terms in Eq. |4.1C| are what is usually defined as the renormal¬ 
ized vertex. The last two terms are one half of the renormalized fermion and 
antifermion self-energy: 


Mt>) = Mt>) - (± - m)5Z 2 - Am 


(4,11) 


where 5Z 2 and Am are coefficients of the counterterms in the Lagrangian den- 
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sity (the ± sign corresponds to the particle/antiparticle case; one half because a 
fermion propagator connects to two vertices). In the above equation, S 2 stands 
for the bare fermion self-energy: 


v / a s f d A k 1 


7m 


(4.12) 


4n J m 2 k 2 (p — k) 2 — m 2 

Upon evaluation, we can write the result for Eq. |4.12| in the form: 

= (± ft — m)E a (p 2 ) + mE b (p 2 ) 

separating it into a scalar and a spinorial component. With these notations, the 
renormalized vertex correction can be written as: 


= W" + i A Z,( P1 ) n + 1 r“ a z 2 (p 2 ) 


with 


A Z 2 (p) = E a (p 2 ) + 


2 | mE b (p 2 ) — A m 


(4.13) 


(4.14) 


± ft — m 

The counterterm coefficient Am is fixed by the on-shell renormalization condition: 


S 2 (i> = m) = 0 Am = mE 6 (m 2 ) 


(4.15) 


Also, in the on-shell limit, 


AZ 2 (p)\^ m = S a (p 2 ) + 


dE b {ft- ft) 

d ft 


<9E, 




d ft 


= 5Z 2 (4.16) 




where for the last equality we have used the renormalization condition 

dMt>) 


d ft 


= 0 




It is convenient to write the contribution of the fermion self-energy diagrams 
in a form similar to that of Eq. |4.9| . Using the resummed top quark propagator: 

i i(± ft + m) 


± ft — m 


p 2 — m 2 


(4.17) 
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we obtain the following result for A Z 2 : 


A Z 2 {p)= [E a (p 2 ) + 2E ir (p 2 )} + 


^ir{p 2 ) 


m 


(± — m) 


(4.18) 


with 


S ir (p 2 ) = m 


E b (p 2 ) - E b (m 2 

p2 _ rfipi 


The term in square brackets in Eq. [4.18| will multiply the Born cross section. The 
term proportional to (± /p — m) is identical to the like terms appearing in the 
expression for the vertex correction Eq. [h9|. Since E ir is the part of the self-energy 
correction which would be infrared divergent on-shell (which means that it has a 
logarithmic resonant behavior 


Zir(p 2 ) ~ log(p 2 - m 2 ) 


in the off-shell case), we keep this term also. 


4.2 Gauge invariance and corrections to partic¬ 
ular subprocesses 

The partial amplitudes appearing in Eq. Q2.5|) can be directly related to Feynman 
diagrams and are straightforward to evaluate. However, as mentioned before, they 
cannot be directly identified with corrections to particular subprocesses. For ex¬ 
ample, the top - antitop production vertex diagram (Figure [A4| a)) contributes to 
the correction to production vertex, as well as to interference between production 
and decay and even to interference between top decay and antitop decay, depend¬ 
ing on when the top quark propagators are closer to the resonances. Moreover, 
the amplitudes in Eq. (|2.5|) are not gauge invariant one by one, although their 


sum is. 
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For purposes related to gauge invariance, and in order to be able to perform 
comparisons with the on-shell computation, it is desirable to decompose the total 
amplitude into gauge invariant corrections to particular subprocesses, and inter¬ 
ference between these. The aim is to rewrite Eq. ( |2.5| ) as: 

M vg = M prod + Mtdec + Muec + M l ^ d _ tdec + M l ^ d _ ldec + M™ de {_ Uec (4.19) 


with each term being gauge invariant by itself. 

To this end, it is necessary to decompose the amplitudes A4 tt , -M#,... into parts 
which contribute solely to corrections to production, decay, or interference. This 
decomposition will be based on the top and antitop propagator structure of the 
matrix element. Following the prescription introduced in section 3.1.1, products 
of propagators which go on-shell in different regions of the phase space can be 
decomposed as follows: 

1 1 - 1 1 1 1 

Dipt) D(pt-k ) D 0 (p t — k) D(pt) D 0 (pt-k ) D(p t — k) 

with 

D(p) = p 2 — fa 2 , D 0 (p — k) = (p — k) 2 — p 2 (4.21) 


In Eq. fP0| the first term on the right hand side is considered as a contribution 
to the production process and the second one a contribution to the decay process. 
Furthermore, it is convenient to write the result in term of products of gauge 
invariant currents (in a manner similar to [|HJ). For example, the A ibt amplitude 
can be written (using the extended soft gluon approximation): 

2 Vb D{pt) 


Mu(ESGA) = G (k) 2v l 

Air J m 2 Dipt A k) D 0 {p b - k) D{p t -k) 


(4.22) 


where Gp U {k) is the gluon propagator in an arbitrary gauge: 

k„k 7J 


G^(k) = 


—i 


k 2 + ie 




k 2 


(4.23) 
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By using the currents : 


• 6 , [i _ 

J tdec 




Dipt 


■ t, p _ 

Jprod 


i M 

D 0 (pt + k) 


D 0 {Pb ~ k) D(p t - k) 

> I tdec 


(4.24) 


D(pt 


D 0 (pt + k ) 


we get : 


a* 


MuiESGA ) = -0 Mo 


d 4 A; 


^(*0 ( jprod - j idee ) jtdec 


■ t, p b, v 


(4.25) 


47 t j m* 

where the hrst term in parentheses contributes to production-decay interference, 
and the second one contributes to decay-decay interference. 

We can similarly define the currents : 


■ t, fl _ 

Jprod 


■ b, p _ 

J tdec 


Apt 


M Dipt) 


Do{p- b + k) D(pt + k) 


(4.26) 


• t, p _ 

J tdec 




Dipt 


Doipt-k) ’ Jtdec Doipt-k) Dipt-k) 

and the amplitudes for the other two interference diagrams can be written like: 

cv 

M tl (ESGA ) = -± 


o 

«s>. | O. 

G^(k) ( jprod-Jtdec ) jtdec 

(4.27) 

. j 

4^ M ° 

[ 0 G^k) 3,t • 

J LI ( 

(4.28) 


Contributions to interference between subprocesses do not come solely from 
the manifestly non-factorizable diagrams. The diagrams in which the gluon con¬ 
tributes to vertex or self-energy corrections (Fig. 2a)) also contain interference 
terms. Since the decomposition into purely vertex (or self-energy) corrections 
and interference corrections is not unique, we shall present onr approach in some 
detail. 

The amplitude for the vertex correction diagram with off-shell particles can 
be written as : 

a. 


mn = £ / 


d A k ^ A(pi,p 2 ) + k a B a (pi,p 2 ) + k a k l3 C a pip ll p 2 ) _ 

2 T II 149 ___7\ /._ , 7. \ 9 9\ W 


m 


{(pi — k ) 2 — m\) (p 2 + k) 2 — ml) 


(4.29) 



CHAPTER 4. VIRTUAL CORRECTIONS 


58 


A corresponding on-shell approximation for this amplitude would be 

l\D - f—r (h\ v + kaB M’P' 2 ) + k a YC a Yp'nP' 2 ) 

47r J in 2 ^ ^ (—2v\k + k 2 ) (2v' 0 k + k 2 ) 


Tu 

(4.30) 

where p\ and p 2 are some on-shell approximations for p\ and p 2 . Now, we can 
dehne the interference contribution through: 

(hr) 12 = (hr)- + (hr)l^ (4.31) 

Note, however, that (<7T)°| is not unique, since p\ , p 2 are not unique; different 
choices for these momenta would yield different results for (hr)°|. The uncertainty 
which arises is, of course, of order p 2 — m 2 , so it can be neglected in the DPA. 
However, it allows us to choose the following definition for (hT)°|: 


ms = 


(y, s 

47T 


d 4 k 


ITT- 


G^(k) 


Ml r (~ 2 ^) 

A>(Pi - k)D 0 (p 2 + k) 


7 


,k a B Q (p u p 2 ) + k a k/ j C a fj(p\,p 2 ) 


7 


(4.32) 


D(pi - k)D(p 2 + k) 

This choice means that the purely vertex correction (factorizable) part of the ver¬ 
tex diagram can be obtained by simply replacing the off-shell C 0 (p\ , p 2 , 0, mi, m 2 ) 
function appearing in the expression for (hT) 12 with the on-shell, infrared diver¬ 
gent function C 0 (pi,P 2 ,R , \JV\-, \flY) 

Conversely, the interference part of the off-shell vertex correction diagram is: 

Mpi,p 2 ) M r (-2 Pi) 


mV = g /g Gpk) 


77T Z 


7 


Dfa - k)D(p 2 + k ) 


Y - 


D 0 (pi - k)D 0 (p 2 + k) 
(4.33) 


with A(pi,p 2 ) = (fa + mi)T(— fa + m 2 ). For the tt production diagram, in DPA 

7 P (fa + rht) r« (- fa + fh t )Y -»• (2pf) r« (- 2 p t -) 

leading to 

n Mu _ _ _ _ 

t , /X 

dec 

(4.34) 


-K? 7 = £ M °J 0 <^(*0 [(-.M) M +iprod(-jtiY + (Mcf) (M<?) 
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Things are different for the decay vertices corrections, since we have doubly 
resonant contributions which are not proportional to the tree level matrix element. 
In the top decay case, the transformation: 

7 M (A + rh b ) £ w + {A + fht) Y -> (2!p£) fi w + [2 p v t + Yi~ A + rh t )\ 


will lead to: 


j<‘ = - 


Oi s 

47T 


d 4 k 


^ G^ik) j&f [HploAMo + Ml 


t, V 


where 

Ml’ 11 = 


-1 


u(b) £ w + Y l r 7 ,z 0 (- fa + m) 4w-v(b) 


(4.35) 


(4.36) 


D(p t ) D{pt) 

In a similar manner, the interference term M'A^ coming from the antitop vertex 
correction diagram can be written in terms of the currents j^ec , —jjrod , and the 
matrix element M\’ v . 

Finally, the last diagrams to be split into on-shell and interference contribution 
are the top, antitop self-energy diagrams. Using the same approach as in the vertex 
case, we define: 


(A Zf 1 = Z { I G^k) 


a 


d A k 


47 T 




7 


a 


d 4 k 


JZ TUT G^ik) 


A + m t 

D{pt - k) 

2 p\ 




7 


2^7 


u 

A 


DoiPt - k) 

2Pt 


+ m t 


Dipt 


47T 


m 


Doipt - k ) A)(Pi - k) 


1) 


(4.37) 


where the quantity in the the curly brackets is the renormalized top self-energy, 
and the quantity on the second line is the on-shell limit of the quantity on the 
first line. This will lead to the following result for the interference contribution 
coming from the top self-energy diagram: 




a 


d A k 


47T 




*’ ** Mr 


Jprod 




t, fi 


■ t , V 
Jtdec 


(4.38) 



CHAPTER 4. VIRTUAL CORRECTIONS 


60 


and a similar one from the antitop self-energy diagram. 

Now we have all the pieces needed to write down the interference terms. The 
final result is: 

d A k 


M intf = 


(Xa 


GW(A0 (j^ od M 0 + »)j” ec - {j P r 0d Mo - Ml’ »)j^ c 


47t J m 2 

+iuJiLcM„] (4.39) 

It is easy in this formula to identify the production-decay or decay-decay interfer¬ 
ence terms. The currents: 

jprod jprod I prod. ' jtdec Jtdec Jtdec > jtdec j Idee I ld.ee (4.40) 

are conserved, and gauge invariant in DPA (that is, = 0, and the partial 
amplitudes in Eq. |4.39| are independent on the gauge paramenter £ in G[k)). 
Therefore, the total interference contribution as well as the interference between 
subprocesses parts are gauge invariant in the approximation used. 


4.3 Computational Approach 

Once a consistent scheme for evaluating the virtual corrections to the top pro¬ 
duction and decay process (|47T|) has been set up (as described in the previous 
sections), the next step is the implementation of this scheme in a Monte Carlo 
program. In this section we give some details about the technical issues arising in 
the design of such a program, and how we choose to solve them. 

There are two types of quantities involved in the evaluation of the NLO am¬ 
plitude: scalar quantities (form factors), which encode the contribution of loops, 
and spinorial quantities, built from Dirac spinors and operators. For example, the 
contribution coming from the tt vertex correction can be written: 

Ms= u(pt) Ai'+(A + "*<) <5r"„(“ faz + m t ) faw-v(Pb) 


(4.41) 
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or, using the decomposition in Eq. 


i=l,2 


fY tY + c A fY t a 


(4.42) 


(definitions for the quantities appearing in the above equation can be found in 
Appendix B). 

Let’s start by discussing the evaluation of the scalar form factors F % . Rather 
than compute analytic expressions for each form factor, we have chosen to evaluate 
them in terms of Passarino-Veltman (PV) functions [Q. This approach has the 
advantage that we have to compute only a few quantities which contain logarithms 
and dilogarithms: the Bq two-point and C o three-point scalar functions (all the rest 
can be written as linear combination of these functions). In turn, for evaluating 
the PV 2 and 3-point scalar functions, we use the FF routines [[hj . 

To compute the amplitudes corresponding to the interference diagrams, we 
need to be able to evaluate the 4-point and 5-point scalar integrals in Eqs. |4.4|, 
EH O There are no published results (or routines) for the general (complex 
masses) 4-point scalar integrals. We have build such routines for the infrared 


finite T > 0 function by using the general methods described in |I5[]. The results 
of these routines have been checked against analytical results in the soft gluon 


approximation published in [[21 


The 5-point scalar function £ 0 has been computed by reduction to 4-point 
functions, following the recipe in [|19[. The resulting infrared divergent 4-point 
functions have been evaluated using the analytic results published in @. 

Some comments on the treatment of the top width are needed here. One way 
of evaluating the scalar form factors in Eq. 4.42| is to compute the gluon integrals 
in the zero top width limit and introduce the finite width only in terms which 
are divergent on-shell (that is, replace m'f with rhf = — im t Tt in terms like 


log {p 2 t 


■m 


! ); see for example |19|). The difference between this result and the 


one 
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obtained by using the complex top mass in all the terms is of order Y t /m t , therefore 
at about 1% level. This would be acceptable if the radiative corrections would be 
small with respect with the tree level result (as is the case for the W production 
process), but in our case it turns out that the one-loop QCD corrections are of 
the same order of magnitude as the tree level result []. Therefore, order % terms 
are important. Since in the case of real gluon radiation the top width appears in 
all terms, for reasons of consistency we need to keep the width in all terms in the 
evaluation of the virtual corrections too. 

The other elements needed in the evaluation of the amplitude ( |4.42|) are the 
spinor sandwiches R. We compute these quantities using spinor techniques, as for 
the real gluon radiation case. Since this part of the computation is quite complex, 
and hence prone to errors, we have two different ways of performing it. In one ap¬ 
proach, we express the R 's in terms of basic spinor products u(pi, Si)u(pj, Sj ); this 
is the more involved case (in terms of the work done by the programmer), geared 
for implementation in a Fortran routine, and which allows fast computation. The 
other approach uses C++ routines which allow the automated evaluation of gen¬ 
eral spinor sandwiches like 

+ rai)(i >2 + m 2 ). ..u(p',s') 

(To this purpose, we have constructed classes that describe < bra | and | ket > 
spinors, and operators of type Ri ±m ; ; in turn, these classes use the basic classes 
- 4-vector, complex number - defined in the Pandora event generator |37]]). This 
method allows easy evaluation of R expressions (again from the programmer’s 
viewpoint) and is much more resistant to programming errors; but the computa¬ 
tion is slower than in the previous method. Therefore, the main use of the results 
5 the reason the total QCD corrections are of order 10 - 20% is because of large cancellations 
between the virtual corrections and soft gluon real corrections. 





CHAPTER 4. VIRTUAL CORRECTIONS 


63 


obtained from the C++ routines is to check the Fortran results. 


4.4 On-shell DPA 

The issue of interference effects in the production and decay of heavy unstable 
particles has been the subject of extensive studies in the past decade. Maybe the 


most important result is a theorem, due to Fadin, Khoze and Martin |38[, which 
states that these interference effects are suppressed. A stronger version of this 


theorem [|39| claims that NLO interference effects cancel in inclusive quantities 
up to terms of order oT/M. Following the methods used in [ [F| , it is possible to 
define a framework for the computation of interference corrections in which the 
total interference contribution to inclusive quantities is zero. 

In this section we shall discuss this alternative approach (which we shall call 
on-shell DPA) to the computation of NLO corrections to the production and decay 
of unstable particles. Results obtained using this approach have been presented 
for the W pair production case at LEP 11 (||l!|, recently, this approach has 


been also applied in the computation of interference (non-factorizable) corrections 
to the top production and decay process at e + e~ as well as at hadron colliders 

ig. 

The relevant features of this approach are two: first, the amplitudes for correc¬ 
tions to subprocesses are computed in the on-shell approximation. For example, 
the correction to the production process can be written in terms of the on-shell 
amplitude: 


M°; rod = xx (e + e- 

A,A' 


tt(g)) M x {t -> bW + ) M x >(t -> bW~) (4.43) 


where A, X' are the spins of the top quarks. The difference between the above 






















CHAPTER 4. VIRTUAL CORRECTIONS 


64 


amplitudes and the ones used in our computation (Eqs. |4.19| , |3T6| ) is due to non- 
doubly resonant terms, therefore acceptable in DPA. 

The other characteristic feature of the on-shell DPA method is that the inter¬ 
ference due to real gluon radiation is computed by using a semianalytic approach. 
This approach rests on the observation that interference is due mainly to gluons of 
energies of order T f ; therefore, we can use the (extended) soft gluon approximation 
in the evaluation of interference terms. 

There are two stages where this approximation comes into play. First, we 
apply it at the matrix element evaluation level. For example, consider interference 
between the diagrams where the gluon is radiated from the bottom quark and from 
the antitop quark. Using the formulas in the Appendix B: 


da r ^(p b ,p w ,...,k) ~ 2Re 


= \M 0 (pb,Pw, ■ • -)| 2 2Re 


4 PiPb 


E Ml 9 (M 


2 — 2 
p t -m t 


sg^* 


(4.44) 


2 kp b + ie ( p t + k ) 2 — ml (pt + k) 2 — ml* 

The second stage is the treatment of the final state phase space. In the soft 
gluon approximation, we can factorize it: dfib,w,...,g = dQb,w,... x dO g , and perform 
the integration over the gluon momenta analytically: 


—Re 

7T 


dv r ib(p'b,Pwi---) = \M 0 (p? b ,p ' w ,...)\ 2 x 

d 3 k IptPb Pt ~ 1 


(4.45) 


27to; 2 kp b + ie (p t + k ) 2 — ml (pi + k ) 2 — ml* 
where p' b ,Pw, ■ ■ ■ are given by a suitable projection of the off-shell momenta 
PbiPWi ■ ■ ■ onto the on-shell phase space (for an example of how this projection 
might be accomplished see |14|). In the above equation we also have p t = p' b + 
p ' w+ , Pi = Pi T p ' w -, and the integral over gluon momenta is allowed to go to 
infinity (since hard gluons contribute nonresonant terms to the result). 
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The quantity on the second line of Eq. [4.45| can be evaluated analytically, 
through methods similar to those used to evaluate the virtual 4-point functions. 
We will not give the results here (they can be found in [Hj, |21|), but there is 
an important comment to make. If we use this procedure to compute the real 
gluon interference, the total interference obtained by adding the virtual diagram 
contribution (Eq. fO|) to the above result and integrating over the top invariant 


mass parameter is zero. The proof of this statement can be found in |39|. This 
cancellation works also for the other interference diagrams; therefore, in this ap¬ 
proach, the contribution of non-factorizable corrections is zero to the total cross 
section. However, this result depends on two things. First, it requires an inclusive 
treatment of real gluon radiation, with phase space integration extending to in¬ 
finity. Second, both the virtual and the real interference terms have to be treated 
in the soft gluon approximation. 

But, is the use of the soft gluon approximation justified in this case? At the 
amplitude level (Eq. |4.44|) , the answer is yes; the relevant gluon energy, being 
of order Tj, is much smaller than the other momenta involved. However, this 
approximation does not seem to be acceptable for the phase space factorization 
stage of the above approach. Here, problems might arise when we try to perform 
the projection of the off-shell momenta onto the on-shell phase space. The reason 
for this is that there is no single way to perform this projection; therefore, in the 
determination of the the on-shell momenta p^p'w, ■ ■ ■ there is an uncertainty of 
the order of the gluon energy, or T t . Now, being close to the top resonances, we 
are in a region of the phase space where the cross section varies greatly over a 
range of energy of order T t (due to the top quark propagators); therefore such an 
uncertainty is not acceptable. 

To illustrate the dependence of the result for real gluon interference on the 
choice of the on-shell momenta p' b ,p' w , ■ ■ ., let’s presume that instead of projecting 
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Pb into p' b , we also take into account the gluon momentum: Pb + k —* p' b (physically, 
this might be justified by the inclusion of the gluon jet in the bottom quark jet). 
Then, Eq. 4.45| becomes: 

d 3 k 4ptPb 1 


d(J lb(Pb,PwT--) = \-M 0 (p b ,p' w ,---)\ — Re 


7T 


2mv 2 kpb + ie {jpi + k) 2 — mf* 

(4.46) 


The result for the above expression is different from the result for Eq. |4.45| , 
and the difference contains doubly resonant terms. Therefore, in the on-shell 
DPA approach, the result for the interference terms depends on how we perform 
the phase space factorization. A discussion of this dependence for the W pair 
production case can be found in |T^ . 


4.5 Results for virtual corrections and the total 
cross section 


In this section, we present some results on the total cross section for the top pro¬ 
duction and decay process at linear colliders. We take into account the virtual cor¬ 
rections as well as contributions coming from real gluon radiation. Furthermore, 
we study the effect of interference (nonfactorizable) terms on invariant top mass 
distributions and we perform comparisons with results previously published j|D 


In obtaining the results presented in this section, the following set of parame¬ 
ters is used: 

m t = 175 GeV, a s = 0.1, T° = 1.55 GeV, T t = 1.42 GeV, 


where T° is the top width at the tree level, while T t includes QCD radiative 


corrections. 
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We start by looking at the total cross section for onr process. Table 4.1 presents 
results for the following quantities: 

• do : cross section for the tree level process fl2.1|) ; computed in the on-shell 
(narrow width) approximation, using the zero-order top width. 

• cr° s : cross section for the NLO process in the on-shell approximation (com¬ 
puted using NLO top width). 

• a{ a,t : the main (factorizable) part of the DPA approximation to the NLO 
process. This quantity contains corrections to production and decay as 
defined in section 4.2. 

• a] nt 'l : the interference (non-factorizable) part of the DPA approximation to 
the NLO process, as defined in section 4.2. 

We present results for three values of collision center-of-mass energies: 360 
GeV, just above the tt production threshold, 500 GeV, the most common value 
used in linear collider studies, and 1 TeV, which can be relevant for higher energy 
machines. Note that at 360 GeV our results are probably not good, being too 
close to the threshold; however, it is interesting to see the magnitude of the 
nonfactorizable corrections at fixed order in this energy range. 

The NLO cross sections contain contributions from the virtual corrections as 
well as from real gluon radiation. We use a technical cut to separate the infrared 
from the real gluons e = 0.1 GeV; the results are independent of the choice of this 
parameter. No physical cuts have been imposed on the final phase space. 

There are several comments to make concerning these results. First, let’s 
compare cr° s with a( act . The difference between these two quantities is due to 
non-doubly-resonant terms, therefore it could be expected to be small. This is 
indeed the case at 500 GeV; but at 1 TeV, this difference is about 6% of the cross 
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Tabic 4.1: Total cross sections for top production at linear colliders (measured in 
picobarns), with no cuts on phase space. 

section. The reason is that in obtaining these results, we have integrated over 
the complete kinematic range available for the top quark invariant mass (that is, 
mb + rriw < \f$ < W — (mb + mw)), so we get contributions from regions of 
the phase space where the top quarks are far off-shell and non-resonant terms are 
important. 

In Tabic 4.2 we present the cross section results obtained with a cut on the t, t 
invariant mass |ypf, ^/p| — m*| < 15 GeV. The difference between the two results 
for the main terms a° s and a{ act is small in this case. Note that, since in the 
on-shell approach p^pj = t?t|, cr° s in Table 4.1 and 4.2 contains a factor which 
simulates the effect of cuts (either from kinematic constraints or imposed ones) 
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360 GeV 

500 GeV 

1000 GeV 

cr° s 

0.682 

0.627 

0.175 

fact 

W 

0.670 

0.629 

0.178 

intf 

-0.034 

-0.007 

-0.002 


Table 4.2: Total cross sections for top production at linear colliders (measured in 
picobarns), with cuts on the top, antitop invariant mass. 

on the top invariant mass. Note also that these cuts are not imposed ad hoc , but 
they arise rather naturally in the process of defining a t,t production event; it 
makes sense to require that the reconstructed mass of the b , W pairs is close to 
the top mass in the definition of such an event. In this context, it is also worth 
noting that the contribution coming from the phase space region where either the 
t or t is far off-shell (more that ten times the width) is quite sizable (around 5% 
of the total cross section for CM energies greater than 500 GeV). 

Another quantity of interest is the differential interference cross section as a 
function of the top invariant mass. Even if the total interference contribution to 
the cross section is small (at about 1% level), it can have larger effects in differ¬ 
ential distributions since it can be positive in certain regions of the phase space 
and negative in others. In particular, it can be important in the reconstruction 
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Figure 4-3: The relative nonfactorizable correction to the invariant mass distri¬ 
bution; the solid line is the contribution of terms proportional to the tree level 
amplitude, while the dashed line contains also the AR terms in Eq. [4.39. 


of the top invariant mass; since da± tends to decrease as Jpf increases, the net 


effect would be to shift the position of the Breit-Wigner peak to smaller invariant 
mass values. This effect can be quantified by the following equation: the shift in 
the mass is 



(4.47) 


where M t = \Jpt, and S n f is the ratio of the non-fa.ctorizable (interference) part 
of the cross section to the Born cross section: 


da\ ntf 



In Figure |4-3| we present the differential distribution for the relative non¬ 
factorizable correction 5 n f(M t ) at center of mass energy 500 GeV. The dashed 
line is the result which takes into account the full interference corrections in Eq. 
4.39| ; the solid line is obtained by taking into account only the terms proportional 
to the Born amplitude. Note that, although the contribution of the M\ terms in 
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Figure 4-4: Real gluon interference: the cr&f term. The solid line corresponds to 
the semianalytic approach; the dashed line is obtained through numerical evalu¬ 
ation with Mt = yjplw'i the dotted line is obtained through numerical evaluation 
with M t given by Eq. |4.48 . 


Eq. |4.39| to the total cross section is very close to zero, they have a sizable effect 
on the differential distribution in Fig. [4-3| . Using Eq. |4.47| , we conclude that 
the shift in the position of the peak in the top invariant mass distribution due to 
interference effects is very small (of order of a few MeV). 

The results in Table 4.1 and 4.2 indicate that the contribution of interference 
terms to the total cross section is of order 1%, in agreement with the T t /mt order 
of magnitude expected from naive arguments. However, it is not zero, as implied 
by results presented in ||0[ , which use the on-shell DPA method. We have argued 
in section 4.4 that this difference is due to the way in which the radiation of real 
gluon with energies of order is treated. In Figure [1 1| . we present the results 
for the real interference between the diagram where the gluon is radiated from 
the bottom quark and the diagram where the gluon is radiated from the antitop 
quark. The solid line is the result of the semianalytical approach described in 
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section 4 . 4 . The other two lines are the result of the exact off-shell computation 
(where the integration over the gluon momenta is performed numerically). The 
two lines correspond to two different ways in which the gluon momentum is treated 
in the reconstruction of the invariant top mass. For the dashed line, the gluon 
momentum is ignored in the top mass reconstruction: M t = \Jp bW . Note that in 
this case, the result is quite close to that of the semianalytical computation, which 
is natural, since the gluon momentum is treated in both cases the same way. 

To obtain the dotted line, we have followed a more realistic approach, in which 
the gluon is included in the top mass reconstruction if it happens to be radiated 
close enough to the top quark: 


M t = 


sjplwg if cos9 tg < vr/3 
\Jpl w otherwise 


(4.48) 


Although the total cross section is the same as for the other exact evaluation case, 
the differential cross section differs by quite a bit. 

The total cross section corresponding to the interference term presented in Fig¬ 
ure Balias the value a b i = 0.121 pb for the semianalytical result, and a bb — 0.124 pb 
for the numerical one. Note that, since the diagram set contributing to this inter¬ 
ference term is not gauge invariant, this result has no physical meaning by itself. 
However, from these numbers we can get some insight concerning the evaluation 
of interference corrections. First, note that the contribution of this single dia¬ 
gram is much bigger (about two orders of magnitude) than the total result for 
the interference terms. This means that there are large cancellations taking place 
between the real and virtual interference contributions. This is quite natural, in 
accordance with the discussion in section 4 . 4 ; however, this also means that small 
uncertainty (order percent) in evaluating one of this contributions (the one com¬ 
ing from the real gluon interference, for example) can lead to large uncertainties 
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Figure 4-5: The relative nonfactorizable correction to the invariant mass distri¬ 
bution; comparison between the semianalytical (dashed line) and the numerical 
(solid line) approach. 

in the evaluation of the total interference contribution. 

In Figure |4-5| we present the comparison between the relative non-factorizable 
corrections computed in the semianalytical approximation (dashed line) and the 
complete off-shell approach (solid line). For the purpose of this comparison, we 
consider only the terms proportional to the Born amplitude in the exact computa¬ 
tion, since only these terms are taken into account in the semianalytical approach. 
Note that the total interference cross section integrates to zero in the latter case, 
and, as discussed above, the complete off-shell distribution contains contributions 
that do not cancel in the total cross section. 



Chapter 5 
Conclusions 


In this thesis, we have discussed in some detail the evaluation of next to leading 
order QCD corrections to the top production and decay process at a linear col¬ 
lider. Since a full computation of the NLO amplitudes contributing to the process 
e + e~ —> b W + b W~ is not feasible, we have employed the double pole approx¬ 
imation. In our case this means taking into account only the diagrams which 
contain two intermediate top quarks. Unlike most of the previous treatments, we 
allow for the two top quarks to be off-shell, and include also the corrections due 
to interference between the top production and decay processes. 

This thesis can be roughly split into two parts: the first one deals with the 
radiation of a real gluon in top production and decay, and the second one deals 
with the computation of virtual corrections and the radiation of soft gluons in 
this process. For the real gluon radiation case, we give a method to split the 
contributing amplitudes into parts which can be thought of as associated with the 
top production or decay subprocesses. This allows for the separation of the cross 
section into parts describing gluon radiation in top production, top decay, and 
interference. We discuss the properties of the gluon radiation, and analyze the 
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impact the gluon has on top mass reconstruction. We pay special attention to the 
analysis of interference terms; since the magnitude of these terms depends on the 
top width, they might provide a way to measure this quantity at energies above 
the threshold. Even if this is not possible, observing experimental evidence for 
interference between the radiation at various stage is an interesting goal by itself. 
Although further studies are needed, our analysis indicates that, while difficult, 
this goal is not a priori unreachable. 

For the virtual gluon corrections case, we first discuss the evaluation of con¬ 
tributing amplitudes in the double pole approximation. A parallel is drawn be¬ 
tween our computation and the DPA evaluation of QED corrections to the W pair 
production and decay process at LEP. While there are many similarities between 
the two computations, there are also some differences; maybe the most important 
one is that there are nonfactorizable (interference) corrections no longer propor¬ 
tional to the Born amplitude in the top quark case. Previous analyses state that 
the interference terms cancel out completely in inclusive quantities (like the total 
cross section). We discuss the evaluation of the real gluon interference terms us¬ 
ing analytic methods and we point out the shortcomings of this approach. The 
total magnitude of nonfactorizable corrections in our computation is found to be 
of order 1% of the cross section. We also present results for the total cross section 
and the relative nonfactorizable correction to the top invariant mass distribution. 
The effect of nonfactorizable corrections on the top mass reconstruction is found 
to be very small. 

The computations in this thesis are valid at collision energies above the top 
-antitop production threshold. Since we are interested in differential distributions 
of final state kinematic variables, the approach used to obtain the results presented 
throughout the paper is that of numerical simulations. The amplitudes are eval¬ 
uated using spinor techniques, and the integration over the final state variables is 
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performed using Monte Carlo techniques. This approach has the added advantage 
that it allows for the inclusion of experimentally relevant selection criteria on the 
final state phase space. 

The calculation presented here is entirely at the parton level. For more real¬ 
istic simulations, it is necessary to take into account initial state related issues, 
like initial state radiation (ISR), beam energy spread and beamstrahlung. The 
hadronization of the final state partons also has to be modeled. In order to address 
these issues, we plan to provide an interface of our code with Pandora (a general 
physics event generator for linear collider studies which includes ISR and beam¬ 
strahlung effects) and Pythia (a Monte Carlo which simulates final state parton 
shower and hadronization). 

For the future, we plan to extend our computations by including all the lowest 
order diagrams contributing to our process. Also, the evaluation of NLO correc¬ 
tions to the singly resonant diagrams in the on-shell approximation should prove 
feasible. Beyond QCD, we plan to take into account clectroweak radiative correc¬ 
tions to the top production and decay process, and maybe SUSY corrections too. 
The framework we use for our computations is flexible enough to build on the 
results presented in this thesis with the final goal of constructing a comprehensive 
Monte Carlo for top related issues at future linear colliders. 
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Appendix A 
Spinor techniques 


The standard technique for computing a cross section used to rely on the evalua¬ 
tion of the square of the amplitude with the help of trace formulas. More recently, 
since the complexity of the processes being analyzed has increased, this approach 
has proven to be too cumbersome. To understand why this is so, note that if the 
total amplitude gets N contributions from N different Feynman diagrams, the 
evaluation of the square amplitude through trace techniques requires the compu¬ 
tation of N 2 terms. Thus, in recent years the emphasis has shifted toward the 
evaluation at the amplitude level, using spinor techniques. 

These techniques have been developed by a number of people over a number 
of years (for some examples, see |E|]). We will use the particular scheme proposed 
by Kleiss and Stirling in |3(J; for the sake of completeness, we review here the 
general features. 
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A.l Spinors describing massive and massless 
fermions 


A spinor describing a massless particle of momentum p and helicity A can be 
constructed with the help of two auxiliary 4-vectors k 0 and k\ and the basic 
spinor U-(k 0 ): 

“^ (p) = vm u - x(k ° ] (A1) 

The k 0 and ki vectors have to satisfy the following conditions : 

k 0 k 0 = 0 , kiki = — 1 , k$ki = 0 (A. 2 ) 

k\ is necessary in order to define the relative complex phase of the positive helicity 
k 0 spinor in relation to the negative one : u + (k 0 ) = fciU_{ko). 

The basic elements from which the total amplitude can be ultimately con¬ 
structed are the products of two massless spinors: 


U \SPl) U \ a (P 2 ) 


These basic spinor products can be expressed in terms of a single complex function 

s(pi,p 2 ): 

u+(pi)u-(p 2 ) = s(pi,p 2 ) (A.3) 

u_(pi)u+(p 2 ) = t(pi,p 2 ) = [s(p 2 ,Pi)}* 


U\{pi)u x (p 2 ) = o 

The exact form of function s(pi,p 2 ) depends on the choice for the auxiliary vectors 
k 0 and k\. Using the values suggested in |30|| : 


(A.4) 


k 0 = (1,1, 0,0), = (0,0,1,0) 
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we have : 

sipuPz) = (p? - ip z i)^pl - (p v 2 - With p(p) = ^2 ^k~ 0 (A.5) 

V{P 1 ) V{P2) 

The spinor which describes a massive fermion of momentum q can be similarly 
defined, with an added degree of complexity: since the number of inner degrees 
of freedom is four for a massive fermion as opposed to two for a massless one, we 
have an extra degree of freedom in our definition: 

1 

u(q , s) = —(A+ fa ± m)u-(p 2 ) (A.6) 

m 

(the minus sign in the parentheses above corresponds to the antiparticle case) 
where p\ and P2 are any two momenta which satisfy the constraints: 

pI=pI = 0, Pi+P2 = q • (A.7) 

What this extra freedom amounts to is the liberty to choose along which axis the 
spin of the fermion points, ft can be shown that the spin vector is given by 

s = {pi-p 2 )/m. (A.8) 


An appropriate choice for the vectors p \, jj 2 will allow decomposition in the helicity 
basis, for example. However, since in our case we are not interested in specific spin 
states for the massive fermions (the spin of the b quarks is not observable, and the 
top quark is off-shell), we have chosen a definition which, while not particularly 
relevant physically, it is computationally convenient: 


m , m 

Pi = I- 7T^ k 0, P‘2 = 


(A.9) 


2 qko 2qko 

Now, since fauik o) = 0 and u(ak 0 ) = \/ol u(ko) (from normalization constraints), 
we get: 

/r 4- m fi — m 

(A.10) 


u(q, ±s) = ^^J^ u^iko), v(q, ±s) = ^—^ u±(k 0 ) 


V 2 Q k o 


V 2 Q k o 
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with 



(note that for the antiparticle spinor, the sign of the spin is reversed, as it should 
be, since in the massless limit, the antiparticle is identical to the particle of op¬ 
posite helicity). 

With these definitions, we get the following expressions for the massive spinor 
products: 


u+(qi)u_(q 2 ) = s(< 7 i, q 2 ), U-(gi)u+(g 2 ) = t(q u q 2 ) (A.ll) 



with the s and t functions the same as in Eqs. [A.3| , [A.5| . If one spinor represents 
an antiparticle, the corresponding mass changes sign in the second line of the 
equation above. 

The total amplitude can be written in term of the basic spinor products dis¬ 
cussed above. Terms like 


u{pi) fafrj ■ ■ .u{p 2 ) 


are evaluated using the completeness relations 


Y u(p)u(p) =ft + m, Y v (p) v (p) = j> — m 
A A 


(A.12) 


while terms of the form: 


[“A, (Pi )'^ u \ 2 M [“As (P3)7m m A 4 (P4)] 


can be evaluated using the Chisholm identity; for massless spinors: 

2 [ u \{P2)u\{Pi) + u_x(pi)u_\(,P2)} if Ai = A 2 = A 

(A.13) 


w A 1 (Pi)7 /i 'WA 2 (P 2 )^ = 


0 


if Ai ^ A 2 
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A.2 Spinors describing massless bosons 

The only massless external boson contributing to our process is the gluon. The 
polarization vector for the gluon can be built with the help of an auxiliary vector 



(A. 14) 


Here, A stands for the helicity of the gluon ( + or - ); we take into account only 
states of physical, transverse polarization. The freedom given by the possibility 
to choose k a corresponds to freedom in the choice of the gluon gauge. 

In relation to the gluon gauge choice, it is worth mentioning here an impor¬ 
tant point concerning numerical instabilities. When the gluon momentum and 
the auxiliary vector for the gluon gauge are almost parallel (kk a ~ 0), the normal¬ 
ization factor in Eq. |A.14| will become very large. If the total amplitude is gauge 
invariant, this is of no concern (besides the fact that there will be large cancella¬ 
tions between different terms in the amplitude). However, if the total amplitude 
contains terms which are not gauge invariant, even if they are usually small, these 
terms will be enhanced. In other words, the cancellations between large terms 
won’t be exact anymore. Obviously this is a problem; the solution is either to use 
a completely gauge invariant amplitude in our computation, or make sure that 
the auxiliary vector k a points in different direction from the gluon momentum. 


A.3 Spinors describing massive bosons 


The treatment of massive boson (W) polarization states actually amounts to let¬ 
ting the W decay into two massless particles (an electron and an antineutrino, 
for example) and integrating over the momenta of these particles. The sum over 
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polarizations can be evaluated as follows: 


~~ [ dQ oTa* v ; aC = uJrA^uJrA (A. 15) 
^ 8irm 2 J 

r'j and r 2 being two lightlike vectors (the momenta of the two massless particles) 
which add up to the momentum of the massive boson, and the integral is over the 
direction of rq. 

In our computation, the polarization vector for the W boson stands for: 


/ w — 


—pjw 
. v/2 


YPl 


—i 


Pw — Mw + iM w T w 


\ —'U)W D , ^ 

u{y) ~^f^PL u ( e 


(A.16) 


where Pl = (1 — 7 5 )/2 is the projector on left state helicity. Treating the W boson 
in the narrow width approximation means that, at the amplitude square level, the 
term coming from the W propagator is: 


7r 


S(Pw - O 


(p 2 w - Ml,) 2 + Ml,T 2 w M W T W 
Therefore, we can replace the quantity in Eq. |A.16| with an effective polarization: 


4w — 


9w 


7T 


YPl U-{v) u_(e + ) 


2 y M\yY \y , 

at the amplitude level, while in the differential cross section: 

d 2 p v d 2 p e < 


(A. 17) 


dO\y — 


(2n) 3 2E l/ (2tt) 3 2 E e + 


d{pw ~ Pu - Pe+) 


Pw= M w 


(A.18) 




















Appendix B 

Amplitudes and cross sections 
formulas 


In this section, we will present some formulas for amplitudes and cross sections. 


B.l Tree level amplitudes 


For the tree level diagrams corresponding to the lowest order process in Figure 
2-l[ the amplitude can be written as : 


A A -( \ L ftt + mt -k + m t 

Mo — u(p b ) fiw+— -— r 7 ,^o —2-“ T~ fcw- v \Pb) 

Pt ~ m t Pi ~ m t 


with the following notations: 


(B.l) 


r 7,^o = [u(P2)(ieY)u{pi)} 


9\u> 

q 2 


u(p 2 ) ( - A’f )) t,(p,) 

cos 9 W * 


(- ieQtrf) + 


9fu> 


—i 


. 9w 


cos 9 


w 


J qi-Ml + iMzTz 
Y\(V‘ - -4V)) 


X 


(B-2) 
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(the first line represents the contribution of the photon exchange diagram, while 
the second line the contribution of the Z 0 boson exchange diagram). Here g w 
is the weak coupling constant, 9w is the Weinberg angle, and Q t = 2/3 is the 
electric charge of the top quark in units of positive electron charge e. p\ and P2 
are the momenta of the initial state electron and positron; q 2 = {p\ + p- 2 ) 2 is the 
square of the total energy available for the process. The couplings of the electron 
and top quark to the gauge boson Z 0 are given by the vector and axial coupling 
parameters: 

V e = ~ + 2sm 2 6 w , A e = —^ 

We shall usually denote by AA the amplitudes multiplied by the denominators of 
the fermion propagators; for example: 

Mo AA Q 9 — 9 9 _ 9 

Pt - m t Pi - m t 

Throughout the paper, we also use the definitions: 

Pt = Pw+ + Pb , Pt = Pw- + Pi , T ~ n t = m t - im tFt , 


and k is the momentum of the gluon. 

For the process with a real gluon radiated, the amplitudes are obtained by 
inserting the corresponding gluon polarization vector in the appropriate place, 
taking into account the extra fermionic propagator, and modifying the momenta 
in the appropriate way. For example, for the diagram with the gluon radiated 
from the top quark: 


,, - ( \ l fct + mt, ■ j %+Jk + m t ~ fa + mt 

M t = U(p b ) C-W+— - —i-Ws fig) 7 - —TV2 - ~~2 


Pt — ys (pt + k ) 2 — rrit 1,Z ° p^ — mfi 


fiw~ v (Pb) ( B - 3 ) 


while for the diagram with the gluon radiated from the antibottom quark: 


M- h 


u(pb) fiw+ 


fa + mt 
2 - 2 
Pt ~ 




~ Pi~ fik + mt 
C Pi + k ) 2 - 


fiw- 


~ fa~ jk + m b 
C Pb + k) 2 - m 2 b 


{-fas fig)v(pl) 

(B.4) 
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Using the decomposition of the top propagator products described in section 
3 . 1 . 1 , we obtain the following expressions for the gauge invariant amplitudes in 


: 


Alprod 


Mltdecay 


Mltdecay 


M t 

2 p t k 

M t 

2 p t k 

Mt 

2 pik 


Mt\ 

1 1 

2 pik ) 

Dipt) Dipt) 

M b N 

\ 1 1 

2 pbk; 

/ D(p t + k) Dipt) 

M~ b \ 

1 1 1 

2p- b k) 

Dipt) D(pt + k) 


(B.5) 


=,2 


m 


where D(p) = p 2 

We also give here the amplitudes in the (extended) soft gluon approximation. 
The evaluation in this case proceeds as in Eqs. |3.8| , |3 .1 0| , and we drop the jt g /c 
terms in the denominators of Eqs. |3.9| , |3. 1 1| . Then 


M { L ilUm . y (ESGA) = (-i 9 .) (2e,p t ) M„ (B.6) 


and, evaluating the other amplitudes in a similar manner, we get: 

Mprod(ESGA) = (-ig s e gil ) M 0 ( ^ ~ ( B - 7 ) 

MtdecayiESGA ) = (-zg s ej Mo D ( pt + k) z^yy 

Mtdecay(ESGA) = {-zg s ej Mo + yyyy) ^yy ^T^y 

Neglecting the gluon momentum in the denominator of the propagators in the 
above equations (a valid approximation if the gluon energy is much smaller than 
the top width) we obtain the expression (Eq. |2T9|) used in the evaluation of the 
infrared singular part of the cross section: 

MZ = (~ig,e m ) Mo (jL - A) (B.8) 
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B.2 Virtual corrections amplitudes 


The amplitude for the general vertex correction in Figure |4- 1| can be written as: 


= 


(X. s 

An 


d A k 1 u A - jk + mi 


in 2 k 2 ^ (pi — k) 2 — m\ ' 1 ’ ( p 2 + k ) 2 — ml 

Upon evaluation of the integral, the result can be written in terms of two sets of 
form-factors; one for the vectorial part of the vertex correction, one for the axial 
part. The number of form factors needed depends on the specific constraints on 
the process; in our case, when the momenta p\ , p 2 are off-shell, we need eight x 2 
form factors: 


r{c v +c A1 b ) 


? 2 — Jk + rn 2 


lu (B.9) 


SK = j- E (CvFYTy* + C A F A T^) (B.10) 

2 = 1,8 

The definition of these form factors depends on the choice of the spinorial 
elements in terms of which the result is written. In our case, we shall define: 







sr'j 

a s 

~ 4 ^ Cv 

X 




(B.ll) 


[ 

Pi 



F i 

+ 



ryl 1 



F? + 

(A - 

- rni) 

Pi 




+ 

(A 

- m 1 ) 

ryV 



F4 + 



Pi 

(- A - 

- rn 2 ) 

F r 

+ 



n ( 

- A ■ 

-m 2 ) 

FI + 

(A - 

- mi) 

Pi 

(- A - 

- m 2 ) 

f r 

+ 

(A 

- m 1 ) 

r ( 

- A ■ 

- m 2 ) 

Fl\ 


for the vectorial part of the vertex correction; for the axial one, replace p±, 7 ^ 
with pi 7 5 , 7 ^ 7 5 in the expression above. This definition has the advantage that 
when the particle i is on-shell, the terms which contain the ± fei + rn t drop out. 
Also, we have made use of the fact that, if hT M multiplies A M in the full matrix 
element can be thought of as the polarization vector of the weak gauge boson 
in diagram [4-l|) , then (77 + p 2 ) /J 'A /t = 0 . 

We shall evaluate and write the results for the form factors in terms of Pas- 
sarino - Veltman functions (for the definition of these see section Appendix C): 
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F\ = 4[ mi(Ci2 — C \i + C 23 — C 21 ) — m 2 (Ci 2 + C 23 ) ] (B.12) 

if = —2[ 2 PlP2 (C 0 + C u ) + 2(C 24 - 1/4) + 25 0 12 - 1 

—mim 2 C , n +Pi(—Cu + C'12) ~ p\ 0 12 } 

if = -4(0, + 2C n - C 12 + C 21 - C 23 ) , F 4 v = 2m 2 (C 0 + Cn) 
if = —4(C 0 + C n + C 23 + C 12 ) , if = 2 , m 1 (C 0 + C n ) 
if = 0 , if = 2(C 0 + C n ) 


For the computation of the axial form factors, we can shift y 5 in Eq. m to 
the right, and then perform the same evaluations as in the vectorial case. In the 
result for the form factors, this amounts to changing the sign of m 2 , m 2 —» —m 2 , 
and multiply the F 5i .. ., F 8 with (—1). Then : 


if = if + 8m 2 (C 12 + C 23 ) , if = if - dmpu^Cn (B.13) 


pA _ rpV pA _ _pV 
r 3 — r 3 1 r A — ^4 



p A — _p^ p A — —F v 

r 7 ~ r 7 i r 8 “ r S 


The result [4.18| for the fermion self-energy corrections: 



can be similarly written in terms of Passarino-Veltman functions: 


Z«(p 2 ) = ^(l + 2B l (j, 2 ,m 2 )) 


(B-15) 
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with 

A B n (p 2 ,m 2 ) = B n (p 2 ,rh 2 ) - B n (rfi 2 ,m 2 ) , n = 0,1 
If we further define the Xq,Xi form factors through: 

ry r ~\ 

AZ 2 (p) =2 -L Xo (p 2 ) + X 1 (p l )(±^-m) 

47T l j 

the renormalized vertex correction in Eq. |4.13| : 

ST^ n = ST“ + i a z 2 ( Pl ) r“ + i r- a z 2 fe) 


(B.16) 


(B.17) 


can be obtained by making the following redefinitions of form-factors in Eq. |B.12 

F 2 v '* + Xo(pl)+Xo(p 2 2 ) 


B.13: 


jpVA . zA ,r A 


(B.18) 


f v c a - Fr + xM) 

f? a - F^ + XM) 

These general results are easily translated for the specific cases which appear 
in onr computation. Take 


n V,A 


F la = FF A (pi =p t , p 2 = Pi, m l =m 2 = m t ) 
for the correction to the top, antitop production vertex, 

F ub = F ^ A (pi = p h , p 2 = -p t , mi = m b , m 2 = m t ) 


nV.A 


for the correction to the top decay vertex, and 

= = ~Ph P 2 = Pi, mi = m, m 2 = m b ) 

for the correction to the antitop decay vertex. Decomposing the top production 
and top decay interaction vertices in a vectorial and an axial part: 

r 7,z» = y‘ - -4' U&ltf (B.19) 
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1 'V,A* 1 „,5 


4w+,w- — 2 £w+,w~Th 2 £w+,w-TtD 
(see also Eqs. |B.2|, |A.17|) we can write: 

■M« = fz £ [cf C,t Uu + Ci F& T& 

i 

= [N f L t L + Ca F* t „ T* a 

i 

A ~A_ ^ s \ ' I "ritb rpV rpV , /~itb t^A rpA 

A4tb ~ 4^ r V li ’ fb + ° A J i,tb 


(B.20) 


where 


T iJi A) = u(p b ) Av+(A + rn t ) r h + mt) 4w-v(Pb) ( B - 21 ) 
T Sb A) =u(p b ) {ih + mt) r 7iZo (- fc+mt) 4w-v(ph) 


rp(V,A) 


i,ib = u (Pb ) 4w+ U>t + mt) r 7 ,z 0 (- fa + m t ) v (Pb) 


and 


/^itt _ t rt _ At /~itb _ (~ttb _ ^ /^tb _ r~itb _ ^ 

La V — V , La A — —/l , Ly — Uy — “ , La A — La A — ~~ . 


B.3 Cross sections and color factors 

The cross section for the generic process (P 1 P 2 ) —> ( PiPj ...) can be written as: 

( B, 2 ) 

where dDf is the differential volume element in the final state phase space: 

d 3 pf 


dn ‘ ,JL. w 32E f 


np> 


(B.23) 


The sum over spins in Eq. [B.22| is performed over the spins of the final state 
particles. If we consider unpolarized electrons in the initial state, we sum over the 
helicities of these particles as well and multiply by a 1/4 average spin factor. 
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In the process under consideration here, the initial collision takes place in the 
center-of-mass frame. Therefore E 1 = E 2 = W/ 2, where W is the total available 
energy for the collision. Also, since the energies involved are much bigger than 
the electron mass (which is considered zero for all purposes in our computation), 
the speed of particles in the initial state is c. 

Actually, in the formulas for the amplitudes presented in the previous sections 
there is an element missing. Since our final state particles include quarks and 
gluons, the amplitude for the process will contain a color index as well. However, 
because the experimental observables are colorless quantities, in the computation 
of the final cross section we have to sum over these color indices. 

Since the initial state is colorless, the bottom and antibottom quarks in the 
final state of the lowest order process have the same color: 

M° -> M° ij = M 0 8 ij 

(here i would be the color index of the bottom quark, and j the color index of 
the antibottom quark). Summing over the colors in the cross section will give us 
a factor of 3: 

= 3 \m°\* 

In the case of a real gluon in the final state, taking the gluon coupling to a pair 
of quarks with color indices i,j to be — ig s T )“• , we have: 

M rg (M r % = M r9 T“- 

where a is the color index of the gluon. The matrices form a representation 
of the gauge group of QCD - SU(3). In order to perform the sum over the color 
indices, we use the relation: 

Tr[T a T b ] = ±5 a , b 
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then: 


a r9 ~ £ \M r9 \ 2 T“ (T^Y = Y \M r9 \ 2 Tr [T a Ta ] = 4 \M r9 \ 2 

a,i,j a 

The cross section in this case acquires a factor of 4 clue to color summation. 

In the case of virtual corrections, all the contributing amplitudes get a factor: 

M vg -> (M V9 )ij = (AT 9 ) Y T k,i T kj = C F M v9 S i:j 

a,k 

Here k is the color index of the quarks inside the gluon loop. For the sum over 
color indices in the equation above, we have used: 

TV 2 — 1 4 

E T T% = C F S tj , Cf = - = 3 for N = 3 

In the cross section, the virtual gluon contribution gets a factor of 3 from the sum 
over the color indices of the quarks in the final state. Note that the C F factor is 
not included in the expressions for AT 1 ’ 9 throughout the paper. 



Appendix C 

Passarino-Veit man functions 


The method of Passarino-Veltman (PV) functions [|33| provides a systematic way 
of evaluating the tensor integrals appearing in the computation of loop corrections. 


For example, the result for the vertex correction integral |B.9| can be written in 
terms of the following tensor functions: 




d A k {1 ,¥\k^k v } 

in 2 (k 2 — m 2 ) ((k + pi) 2 — m 2 ) ((k + pi + p 2 ) 2 — ml) 


Following the Passarino-Veltman method, these functions can be written in term 
of scalar functions: 


= P tC u +P^Ci2 (C.2) 

C' W = P 1 P 1 C 21 + P 2 P 2 C 22 + {P 1 P 2 + P 2 P 1 ) ^23 + 9^ U C24 

Furthermore, the C tJ functions in the above equations can be written in terms of 
the scalar one-, two- and three-point integrals A°,B° and C°. For completeness, 
we give here the formulas used to compute these functions (these formulas can be 
easily derived by multiplying the Eqs. |C.2| by p i M , p 2 tl , g, JV and solving the resulting 
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systems of equations; they can also be found in p3|| ): 

j(23) + /i C° 



= X 


f _ ^ I 1 K?( 23) 

624 “ 4 + I^° 


-l 


( w(13) v2\ 

a 0 a 0 


\ Bo 


( 12 ) 


B, 


(13) 


1 


+ he 0 ) 


+ i™;c° - i(/,C„ + hCn) 


f C 21 ) 

= A- 1 

V C23 ) 

l 

1 C 23 N 

= A -1 

V C22 ) 

l 


( r( 13 ) 


B'r + B { i 3) + f\C n - 2C 24 

b (12) _ B (VS) + /aCu 
^(13) _ S (23) + /iCi2 

^ -B? 3) + hCi2 ~ 2C 24 


Here we have used the following notations: 

/ 


X = 2 


Pi P 1 P 2 
\ P 1 P 2 p\ 


fi = ml-ml- pi 

f 2 = m \ -ml - (p 1 + p-2 ) 2 + p 2 


H {0 ’ m} (pi, mi,m 2 ) 
B^(pi, mi,m 2 ) 


cRk 


{1 


J m 2 (k 2 — ml) ((k + pi) 2 — m^) 

= PiBi(pl,ml,ml) 


2plB 1 (pl,ml,ml) = 


A 0 (ml) - A 0 (ml) + fiB 0 (pl, ml, ml) 


(C.3) 

(C.4) 

(C.5) 


(C.6) 

(C.7) 

(C.8) 


^i 12) = B n (pi,m?,m|) (C.9) 

B { n 3) = B n ((p 1 +p 2 ) 2 ,ml,ml) 

B {23) = B n (pl,ml,ml) , n = 0,1 



and, finally, 


(C.10) 
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The ultraviolet divergent integrals in the above expressions are evaluated with the 
help of the dimensional regularization method. 

The above equations follow the standard definition of the PV functions; we 
shall denote the C functions in Eq. |C.1| , |C.2| by m\, m 2 , m 3 ) and 

the C functions in Eqs. |C.3| , |C.5| , |C.5| by The C function 

used in Appendix B are defined by: 

{1 ,k^,k^k v ) 


C {0,^} = fd^_ 


in 2 ( k 2 + ie) ((k — pi) 2 — m\) ((k + p 2 ) 2 — ml) 


(C.ll) 


CT = —p\Cu + (p 2 — P1YC12 (C.12) 

C^ U = P 1 P 1 C 21 + (p 2 - PiT(P 2 ~ Pl) U C 22 - [Pi(P2 - PlY + (P 2 ~ PlYPl] C 23 


therefore: 


= C^’^(-pi,p 2 -pi,0,mi,m 2 ) 
Co, Cij = C°, Cij(pl, (pa - Pi) 2 , 0, ml, ml). 


Moreover, the function B q 2 in Eq. |B.12| is given by: 


B^ 2 = B"((p 1+ p 2 y,ml,ml) 


and in Eqs. |B.15 


(C.13) 


(C.14) 


B n (p 2 , m 2 ) = B n (p 2 , 0, m 2 ) , for n = 0,1. 


(C.15) 



















